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ABSTRACT 



Aims. Cosmic shear is a powerful method to constrain cosmology, provided that any systematic effects are under control. The intrinsic 
alignment of galaxies is expected to severely bias parameter estimates if not taken into account. We explore the potential of a joint 
analysis of tomographic galaxy ellipticity, galaxy number density, and ellipticity-number density cross-correlations to simultaneously 
constrain cosmology and self-calibrate unknown intrinsic alignment and galaxy bias contributions. 

Methods. We treat intrinsic alignments and galaxy biasing as free functions of scale and redshift and marginalise over the resulting 
parameter sets. Constraints on cosmology are calculated by combining the likelihoods from all two-point correlations between galaxy 
ellipticity and galaxy number density. The information required for these calculations is already available in a standard cosmic shear 
data set. We include contributions to these functions from cosmic shear, intrinsic alignments, galaxy clustering and magnification 
effects. 

Results. In a Fisher matrix analysis we compare our constraints with those from cosmic shear alone in the absence of intrinsic 
alignments. For a potential future large area survey, such as Euclid, the extra information from the additional correlation functions can 
make up for the additional free parameters in the intrinsic alignment and galaxy bias terms, depending on the flexibility in the models. 
For example, the Dark Energy Task Force figure of merit is recovered even when more than 100 free parameters are marginalised 
over. We find that the redshift quality requirements are similar to those calculated in the absence of intrinsic alignments. 
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1. Introduction 



On its way to Earth the light from distant galaxies is contin- 
uously deflected by the matter density inhomogeneities which 
it passes by. This induces distortions of the shapes of the 
projected galaxy images on the sky, causing modifications of 
these shapes of order 1 %, known as weak gravitational lens- 
ing, or cosmic shear. Hence, to measure this distortion effect 
on distant galaxies, one requir es st atistical me t hods; see e.g. 
Bart elmann & Schneider! (1200 ll) and ISchneiden d2006l) for de- 
tailed reviews. 



Detected in 2000 dBacon et al.1 120001: iKaiser et al.1 l2000t 
van Waerbeke et aD l2000t IWittman et all 2000), cosmic shear 
has since rapidly evolved into a mature technique that pro 
duces increasingly stringent constraint s on cosmologi c al pa 
ra meters (see e.g. Jarvis et al.1 20061 Hoekstra et al. l200< 



ii 



Semboloni et al. 2006; Hetterscheidt et al. 2007; Beniamin et al 
20071: iMassev et al.l 120071: ISchrabback et al.l 120071: iFu et ail 
2008; Schrabb ack"etal.ll2010h Probing both the evolution of 



STARRS0, DEfl LSSlQ JDEIvfl and Euclid^ will produce 
deep imaging over large fractions of the sky and thereby yield 
unprecedented insight into the properties of dark matter, dark 
energy and gravitation (e.g.lTakada & J ain 2004; Refregier et al. 
l2004ll2006LlKitching et aDl2008aHThomas et aljl2009h . 

The great statistical power of cosmic shear demands a care- 
ful assessment of possible systematic errors that might bias the 
results. A serious limitation may arise from a physical systematic 
caused by the intrinsic alignment of galaxies. The matter struc- 
ture around galaxies can modify their intrinsic shape and their 
orientation. Firstly, this can result in correlations between the 
intrinsic shapes of galaxies which are close both on the sky and 
in redshift (intrinsic ellipticity correlations, or II correlations). 
Moreover, a dark matter halo can intrinsically align a physically 
close galaxy in the foreground and at the same time contribute to 
the lensing signal of a background object, which i nduces gravita- 
tional shear-intrinsic ellipticity correlation (GI, Hirata & Seliak 
12004 . 

Detailed investigati ons have been pe r formed on the align- 
ment between haloes dCroft & Metzlerl l2000t iHeavens et all 



structure and the geometry of the Universe at low redshifts, it is 
complement ary to observations of the cosmic microw ave back- 
ground (e.g. iTereno et al.ll2005l iDas & Sper gel 2009) and con- 
sidered as potentially the most powerful method to pin down 
the properties of da rk energy dHul 120021; lAlbrecht etani20Q6l 
iPeacock et al.l l2006h . Upcoming and future surveys like Pan- 



1 Panoramic Survey Telescope & Rapid Response System, 
http: //pan- Starrs . ifa.hawaii . edu 

2 Dark Energy Survey, https : / /www . darkenergysurvey . org 
Large Synoptic Survey Telescope, http: //www. lsst . org 

4 Joint Dark Energy Mission, http : // j dem . gdf c . nasa . gov 

5 http : //sci . esa . int/ science-e/www/area/ 
index . cfm?f areaid=102 
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2000 
2001 



Lee & Peril 120001: . . 
Jing 2002; Mackev et al. 2002 



Catelan etaT] 12001: Crittenden et al| iJohnston efail 120071: 



Hirata & Seliak 2004; 
Bridle &Abdallall2007t ISchneider & Bridldl2010l) . as well as 
the alignment of the spin or the shape of a galaxy with its own 



tne alignment of tne spin or tne snape of a galaxy witn its own 
dark matter halo (e.g.lPen et al.l2000f van den Bosch et all2002t 



Oku mura et alj 120091; lOkumura & Jind 120091: IBrainerd et al.1 
2009; see also Schafer 2009). Intrinsic alignments have 



also been investigated observationally, where non-vanishing 
II and GI signals have been detected in several surveys 



dBrown et all 120021 iHevmans et all 20041 


Mandelbaum et al.l 


2006;lHirata et alj|2007l IBrainerd et al.ll2009 


). 



The results of both theoretical studies and observations show 
large variations, but most are consistent with a contamination of 
the order 10 % by both II and GI correlations for future surveys 
that further divide the galaxy sample into redshift slices (cosmic 
shear tomography). Hence, the control of intrinsic alignments in 
cosmic shear studies is crucial to obtain unbiased results on cos- 
mological parameters. Accurate models would solve the prob- 
lem, but progress is hampered due to the dependence of intrin- 
sic alignments on the intricacies of galaxy formation and evolu- 
tion within their dark matter environment. Currently, the level of 
models i s crude, and partly only m otivated phenomenologically 
(see e.g. ISchneider & Brid le 2010 for recent progress). 

The II contamination can be controlled relatively eas- 
ily by excluding close pairs of galaxies from the analy- 
sis dRing & Schneiderll2002l 120031; IHevmans & Heayensll 2003: 
iTakada & Whitdl2004l)JJoachimi & Schneider! d2008ll2009l) in- 
troduced a nulling technique which transforms the cosmic shear 
data vector and discards all entries of the transformed data set 
that are potentially contaminated by the GI signal. While this 
approach only relies on the well-kn own redshift dependence 
of gravitational lensing, I Kind (120051) projects out the GI term 
by making use of template func tions. Furthe r more the work 
by iMandelbau m et all ([2006) and iHirata et afl ([2007) suggests 
that the intrinsic alignment may be dominated by luminous red 
galaxies which could be eliminated from the cosmic shear cat- 
alogues. All these removal techniques require excellent redshift 
information, and still they can cause a significant reduction in 
the constraints on cosmology. 

Deep imaging surveys not only provide information about 
the shape of galaxies, but allow in addition for a measurement of 
galaxy number densities, as well as cross-correlations between 
shape and number density information. This substantial exten- 
sion of the set of observables increases the cosmological infor- 
mation to be extracted and, more im portantly, enables one to in- 
ternal ly calibrate systematic effects dHu & J ain 2004; Bernstein 
2009). By adding galaxy number density information one adds 
signals that are capable of pinning down the functional form of 
intrinsic alignments, but one also introduces as another system- 
atic, the galaxy bias, which quantifies the lack of knowledge 
about how galaxies, i.e. the visible baryonic matter, follow the 
underlying dark matter distribution. 

It is the scope of this work to elucidate the perfor- 
mance of a joint analysis of galaxy shape and number den- 
sity information as regards the ability to constrain cosmo- 
logical parameter in presence of general and flexible mod- 
els of intrinsic alignments and galaxy bias. In doing so we 
incorporate several cosmological signals which have been 
considered before as promising probes of cosmology them- 
selves, includ ing galaxy clustering from photometric red- 
shift s urveys dBlake & Bridld 120051: iDolnev et al J 120061: IZhanl 
2006; iBlake et all 120071: iPa dmanab han et all 120071) galaxy- 
galaxy lens i ng (e.g. [Sc hneider & Ria 1 1 997t [Guzik & Seliak 
1200 1L I2002t ISehakl l2002t ISeliak et al.1 120051: lYoo et alJ 12006: 



_ iCacciato et al 20091) and lensing mag- 
nification dBroadhurst et alJ 11995b IZhang & Pen! 12005112 006: 
van Waerbekd2 010). We follow the ansat z outlined in | Bern stein 
(20091) and extend the investigation bv iBridle & Kind d2007n 
who considered the residual information content in galaxy shape 
correlations after marginalising over the parameters of two log- 
linear grid models representing the II and GI terms. We quantify 
the cross-calibration properties of the joint set of observables 
and determine the requirements on cosmological surveys to effi- 
ciently apply this joint approach. 

This paper is organised as follows: In Sect. [2] we give an 
overview on the two-point correlations that form part of the 
galaxy shape and number density observables, and we derive 
their explicit form. Two appendices provide further details. 
Section [3] demonstrates how we model the different signals and 
their dependence on cosmology. We introduce a general grid 
parametrisation for the intrinsic alignments and the galaxy bias. 
Furthermore we summarise our Fisher matrix formalism and the 
figures of merit we employ. In Sect.[4]we present our results on 
the dependence of the parameter constraints on the freedom in 
the model of intrinsic alignments and galaxy bias, the character- 
istics of the redshift distributions, and the priors on the different 
sets of nuisance parameters. Finally, in Sect.[5]we summarise our 
findings and conclude. 



2. Two-point correlations from cosmological 
surveys 

Cosmological imaging surveys observe the angular positions and 
the projected shapes of huge numbers of galaxies over increas- 
ingly large areas on the sky. In addition, by means of multi- 
colour photometry, it is possible to perform a tomographic anal- 
ysis, i.e. obtain coarse information about the line-of-sight dimen- 
sion in terms of photometric redshifts (photo-z). From the galaxy 
shapes in a given region of space, one can infer the ellipticity 



7%W + yf(0) + 6® (0) 



(1) 



where the superscript in parentheses assigns a photo-z bin i. 
The observed ellipticity e has contributions from the gravita- 
tional shear jc and an intrinsic shear y\, which is caused by 
the alignment of a galaxy in its surrounding gravitational field. 
Moreover, e is assumed to have an uncorrelated component e rnc j, 
which accounts for the purely random part of the intrinsic orien- 
tations and shapes of galaxies. Note that dTj is only valid if th e 
gravitational shear is weak, see e.g. ISeitz & Schneiderl dl997l) : 
iBartelmann & Schneiderl (12001 1) and for certain definitions of el- 
lipticity. 

Likewise, the positions of galaxies can be used to construct 
an estimate of the number density contrast 



(0, 



,(0 



(2) 



which is determined by the intrinsic number density contrast 
of galaxies n g and the alteration of galaxy counts due to lens- 
ing magnification n m . An uncorrelated shot noise contribution is 
added via n rn( j. In contrast to e (,) (0) the number density contrast 
n w (0) can obviously not be estimated from individual galaxies. 
One can understand n w (0) as the ensemble average over a hypo- 
thetical, Poisson-distributed random field of which the observed 
galaxy distribution is one particular representation. The formal 
relation between the projected number density contrast as used 
in ([2]i and the three-dimensional galaxy number density fluctua- 
tions will be provided below, see (fT2l . 
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As was already noted in Bernstein (2009), (0 and (0 are 
symmetric in the sense that they both contain an intrinsic contri- 
bution and a term caused by gravitational lensing effects. Under 
usual circumstances the correlated part of the ellipticity is dom- 
inated by the gravitational shear, whereas the largest term in (0 
is due to the intrinsic number density contrast. 

Both ellipticity and number over-density vanish if averaged 
over sufficiently large scales. Thus, one considers to lowest order 
two-point statistics of these quantities. Since all real-spac e two- 
point measures are related to the power spectrum (see e.g. lKaiserl 
119921) . we can work in terms of power spectra without loss of 
generality, which is desirable in particular due to a simpler struc- 
ture of the signal covariances in Fourier space. Denoting the 
Fourier transform by a tilde, the power spectrum C^{€) between 
redshift bins i and j can then be defined by 

(xf{t) ^V)) = (2tt) 2 5%\t - t') C^'V) , (3) 

where 5^ is the two-dimensional Dirac delta-distribution, and 
where t denotes the angular frequency, the Fourier variable on 
the sky. The measures x a and x\> can correspond to any of the 
set {yc, Ji, n s , n m }. The random contributions in (0 and (0 are 
not correlated with any of the other measures and only yield a 
contribution to the noise, see Sect. 13.41 

Inserting and into (0, one obtains the complete set 
of tomographic two-point observables which are available from 
shape and number density information 



Cnn (0 — ^gg CO ^gmCO ^gmCO ^mm(0 



■ C»{>({) + c 



mG 



(0 + C™(£) 



(4) 
(5) 
(6) 



see iBernsteinl (120091) . We name signals stemming from galaxy 
shape information by capital letters ('G' for gravitational shear, 
T for intrinsic shear) and signals related to galaxy number den- 
sities by small letters ('g' for intrinsic number density fluctua- 
tions, 'm' for lensing magnification). An overview of the nomen- 
clature of the correlations in (0 to (0 is provided in Table 
Note that (0 and are symmetric with respect to their photo- 
z bin arguments. Hence, if N z bi n denotes the number of avail- 
able photo-z bins, one has iV z binCN z bin + l)/2 observables for ev- 
ery considered angular frequency. In contrast, one can exploit 
N 2 bin ellipticity-number density cross-correlation power spectra 
per {. 

The set of observables in (0 is the one that cosmic shear 
analyses are based on. The shear correlation signal (GG) is 
a clean probe of the underlying matter power spectrum and 
is thus po werful in constraining cosmological parameters (e.g. 
iHul IT999). However, shape measurements incorporate further 
terms stemming from correlations of intrinsic ellipticities (II) 
and shear-intrinsic cross-correlations (IG, or equivalently GI) 
whose contribution can be s ubstantial, but is to date poorly 
known dHirata & Seriakl l2004h . These terms exist because the 
shapes and orientations of galaxies are influenced via the tidal 
gravitational fields of the matter structures in their surrounding, 
which firstly induce correlations between neighbouring galax- 
ies, and secondly cause correlations by determining the intrinsic 
shape of a foreground object and adding to the shear signal of a 
background galaxy. 

Intrinsic galaxy clustering (gg) add a strong signal to the cor- 
relations of galaxy number densities (0, but its use to obtain cos- 
mological parameter estimates is limited due by poor knowledge 



of the galaxy bias (e.g. Lahav & Suto 2004). Gravitational lens- 
ing modifies the flux of objects and thus reduces or increases 
number counts of galaxies above a certain limiting magnitude, 
depending on the form of the galaxy luminosity function close 
to the limiting magnitude. This produces magnification corre- 
lations (mm) and intrinsic number density-magnification cross- 
correlations (gm). The gm correlations occur when a foreground 
mass overdensity (underdensity) contains an overdensity (un- 
derdensity) of galaxies and (de)-magnifies background objects 
along the same line of sight causing an apparent over- or under- 
density of galaxies at higher redshift. 

Cross-correlations between galaxy number densities and el- 
lipticities (0 contain contributions from cross terms between 
intrinsic clustering and shear (gG), intrinsic clustering and in- 
trinsic shear (gl), magnification and shear (mG), and magnifica- 
tion and intrinsic shear (ml). For instance, one expects to find 
gl and gG signals when a mass structure leads to an overdensity 
in the local galaxy distribution and influences the intrinsic shape 
of galaxies at the same redshift or contributes to the shear of 
background objects. The latter is the usual galaxy-galaxy lens- 
ing signal. Because a foreground overdensity can in addition en- 
hance galaxy counts due to lensing magnification, the mG and 
ml signals will also be non-vanishing. The form of all these cor- 
relations will be further discussed in Sect.0 

All non-random terms in and 0, given for a photomet- 
ric redshift bin i, can be related to a source term S, which is a 
function of spatial coordinates, i.e. 



v(0 



(0) 



-r 

Jo 



(7) 



where we defined a weight function w that depends on the photo- 
z bin i (for a similar approach see iHu & Jainll2004l) . Here, x 
denotes comoving distance, and fvXx) is the comoving angular 
diameter distance. If (0 holds for two quantities and x9 , 
their projected power spectrum is given by the line-of-sight in- 
tegral of the three-dimensional sour ce power spec trum Ps a s b via 
Limber's equation in Fourier space (Kaiser| |1992l) . 



r 

Jo 



Mx) 



• X 



(8) 



By identifying weights and source terms for gravitational and 
intrinsic shear, as well as intrinsic clustering and magnification, 
we can derive Limber equations for all power spectra entering 
©-©. 

To compute the equivalent of (0 for the cosmic shear case, 
we first note that in Fourier space the shear and the convergence 
are related by the simple equation kg(() = 7g(() e~ 2m , where <pt 
is the polar angle of t. As a consequence, the power spectra of 
shear and convergence are identical. Therefore, we can equiva- 
lently use the convergence k ( '\0) as the cosmic shear observable. 
It is related to the three-dimensional matter density contrast 6 via 



where the weight is given by 



q W (x) 



3/^Qn 

2 c 2 



I 

Jy 



Ax' P W (X') 



(9) 



(10) 



a(x) J x MX') 

see iBartelmann & Schneider! (120011) : ISchneiderl d2006l) for de- 
tails. Here a denotes the scale factor and p^'Hx) tne comoving 
distance probability distribution of those galaxies in bin i for 
which shape information is available. 
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Analogously to the lensing case, one can define a conver- 
gence of the intrinsic shear field iq (0), which is directly related 
to the intrinsic shear via k\{1) = e~ 2llft . This intrinsic con- 
vergence is a projection of the three-dimensional intrinsic shear 
field Hi, which can be written as 

Kf{0)= dxp & OtiKi(fK(x)0,x) . (ID 
Jo 

see e.g. Eirata & Seljak (2004) for the analogous expression 
in terms of intrinsic shear. Here we have assumed that the in- 
trinsic shear field is - like the gravitational shear field - curl- 
free to good approximation. This holds for instance for t he lin- 
ear alignment model developed in lHirata & Seliakl (120041) . Then 
Ki corresponds to the Fo urier transform of y^ik) as defined in 
ISchneider & Bridld (l2010h . 

Likewise, angular galaxy number density fluctuations n\ -(0) 
are given by the line-of-sight projecti on of three-di mensional 
number density fluctuations 6 g as (e.g. lHu & Jain|l2004T) 

rXhor 

4W= dxp M (x)6 g (MxW,x) • (12) 
Jo 

It is important to note that p®(x) is the same as in (fTOb and ( fTTI ), 
i.e. the number counts are restricted to those galaxies with shape 
measurements, which require a higher signal-to-noise than the 
position determination. In principle, number density informa- 
tion could be obtained for a larger number of galaxies, in par- 
ticular fainter ones. But, to determine the contribution to num- 
ber density correlations by magnification, it is necessary to mea- 
sure the slope of the luminosity function a w at the faint end of 
the used galaxy distribution. We will detail the exact definition 
and the determination of a (i) in Sect. 13.21 Since it is desirable 
to extract the values of the or® internally from the survey, one 
needs to be able to measure fluxes down to values slightly below 
the magnitude limit of the galaxies included in p®(x). Moreover 
galaxy number density measurements may require photometric 
redshifts which are of the same or better quality than for cos- 
mic shear studies, limiting the number of faint usable galaxies. 
Hence, we argue that the choice of identical distance probabil- 
ity distributions for both shape and number density signals is a 
fair assumption. We add the warning that one may have to ac- 
count for selection biases, for instance if one investigates cos- 
mic s hear and magnificat i on effects with the same galaxy sample 
(e.g. lSchmidt et aUfcOOgtlkrause & Hiratall2009l) . 

We write the contribution of magnification effects to the 
number density measurement as 

n%(ff) = 2(a w -l) dxq i0 (x)5(Mx)e,x) (13) 
Jo 

= 2(a®-l)*®(0). 

As several assumptions enter this equation, in particular con- 
cerning the treatment of the slope a (,) , we provide a detailed 
derivation of ( fT3l in Appendix [A] Comparing the projection 
equations (0, (fTTI) . (fT2l i. and ( fT3l to the general form 0, one 
can derive all possible cross- and auto-power spectra in the form 
of the general Limber equation ([§). For completeness we have 
given these Limber equations in AppendixlBl 

Good models of the three-dimensional source power spec- 
tra in the Limber equations (see also the right hand column of 
Table [TJ are unknown except for the non-linear theory matter 
power spectrum Pgg. The distribution of galaxies is expected to 
follow the distribution of dark matter, so that the galaxy cluster- 
ing power spectra should be related to Pg$. However, to date it 



Table 1. Overview on the two-point correlations considered in 
this work. 



measured correlation 


OT~\ TOO 


QT~\ TOC 


shear 


Cqg 


P 66 


intrinsic-shear 


Cig 


Pei 


intrinsic 




Pu 


galaxy clustering 


c gg 


p 

r-gg 


clustering-magnification 


r 

<-gm 


Pg6 


magnification 




Pm 


clustering-shear 


C g G 


Pg6 


clustering-intrinsic 


Cgi 


P* 


magnification-shear 


C m G 


Pm 


magnification-intrinsic 


C m I 


Pm 


galaxy ellipticity (observable) 


c ee 




galaxy number density (observable) 


c 

III! 




number density-ellipticity (observable) 


r 





Notes. Listed are the symbols used for the two-dimensional projected 
power spectra and the underlying three-dimensional power spectra. 



is unknown how much the galaxy clustering deviates from dark 
matter clustering, in particular on small scales. This is usually 
expressed in terms of the galaxy bias b g , which is a function 
of both angular scale k and redshift or line-of-sight distance x- 
Hence, one can write 

P gg (k,x) = b\(k, X ) P ss (k, X ) (14) 
PgsiKx) = b g (k,x) r g (k, X ) Pss(k, X ) , 

where to describe the cross-correlation between matter and 
galaxy clustering, we introduced a correlation coefficient r % in 
the second equality. 

The intrinsic alignment power spectra depend on the intrica- 
cies of galaxy formation and evolution within their dark matter 
environment. Again, precise models of the intrinsic alignment 
have to rely on baryonic physics and are currently not avail- 
able. For symmetry reasons we parametrise our lack of knowl- 
edge about the intrinsic alignment power spectra similarly to the 
galaxy bias as 

Pn(k, X ) = bj(k,x) Pss(k,x) (15) 
Psi(k,x) = bi(k,x) n(k,x) P ss (k,x) 

with the in trinsic alignment bias b\ and correlation coefficient r\ 
(following lBernsteinll2009t) . Although the power spectrum P g \ 
could in principle contain a third, independent correlation coef- 
ficient, we assume that it is sufficient to write 

Pgi(k,x) = h(k, X ) n(k,x) h(k,x) r g (k,x) P 6S (k, X ) , (16) 

i.e. we hypothesise that correlations between intrinsic number 
density fluctuations and intrinsic alignments can entirely be 
traced back to the effects of the intrinsic alignment bias and the 
galaxy bias. This is a strong assumption since instead of intro- 
ducing a fifth completely unconstrained bias term, ( fL6T > estab- 
lishes a link between the galaxy bias and intrinsic alignment bi- 
ases. 

Our equation ( fT6|) is effect i vely in cluded within the last term 
in curly brackets of Bernstein (20 091), eq. (19), and w e have ef- 
fective ly set s gK = in iBernsteinl d2009l) . eq. (35). [Bernstein 
(l2009h fixes his s gK to be a single unknown scalar across the sur- 
vey, stating that we expect this type of cross-correlation to have 
a minimal effect on cosmological constraints. It would be inter- 
esting to check this by comparing results in which s gK is allowed 
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to vary with those in which s gK = 0. However this is beyond the 
scope of this paper. 

We note that the unknown quantity in question is the cross- C^(() - 
correlation between the in trinsic alignment field and the galaxy 
position field (r gK in the iBernsteirj 120091 notation), which is C&&(£) = 
precisely the corre lation measured in the observational con- 
straints papers (e.g . Mandelbaum et al. 2006t lHirata et al.l l2007: 
Mandelba um et al.ll2010l) . Our ansatz states that this is simply 
related to the cross-correlation between the intrinsic alignment 
field and the mass, and the cross-correlation between the gal axy 
position field and the mass (r gK = r 8 r K in the lBernsteinll2009l no- 
tation). These latter two quantities are much harder to measure 
because we do not know the mass field well. Indeed gravitational 
lensing and the observations we discuss in this paper are likely 
the best probe of these quantities. Ultimately these data could be 
used to constrain this additional freedom. 

To shed further light on this question, we now discuss a 
simple toy model to illustrate the nature of our assumption. In 
the simplest case the galaxy distribution may trace the mass 
distribution, and the intrinsic shear field is the gradient of the 
mass distribution. In a toy universe we could displace the posi- 
tion mass field on the sky relative to the other two fields which 
would render the correlation coefficient rj (or r" in the notation 



Xbl 



2(a«-l)C^V) 



fcfxY 



r Xh - p {i) (x) fHx) 

Jo ' " ' 



X 



Jo 
x b 



(21) 
(22) 

(23) 



I 



81 /k(*)"7 b Uk(x) 

d X 72 b 



,X\ Pss 

I 

MxY 



MxY 



x b\ 



Mx) 



<x\ n 



MxY 



,x\ p& 



x h 



MxY 



= 2(a (/) -l)C™ 



"GG V 





1 




t \ 


\fK(x)' X ) 


■*) 


(24) 




(25) 




(26) 



of Bernstein 2009) zero and also the correlation r„ 



(for 



large displacements), leaving a strong correlation between the 
intrinsic shear f ield and the galaxy position field (unity r gK in the 
lBernsteinir2009i notation, thus s gK =1). 

Similar more physical arguments could be made by consid- 
ering stochasticity instead of displacements. However, we might 
expect the physical origin of both the galaxy position field and 
the intrinsic shear field to both lie in the mass field, since most 
of the interactions are mediated by gravity. Therefore it seems 
reasonable to expect there to be no additional interplay between 
the intrinsic shear field and galaxy position fields t hat would 
produ ce further cross-correlations (non-zero s gK in the lBernsteinl 
120091 notation). 

We will not limit the values of the correlation coefficients to 
the interval [-1; +1]. It is formally possible that |r| > 1 if our as- 
sumption about the statistics of the galaxy distributio n, usually 
taken to be Poissonian, is incorrect dBernsteinl 12009) . Treating 
the correlation coefficients as completely free parameters, our 
choice of parametrisation in (TBI) and (15[ is equivalent to mod- 
elling P gg and P g g, or likewise Pn and Psi, independently. 

We insert the parametrisations ( fT4b . dT51 >. and (fT6l i into the 
set of Limber equations and can this way relate all power spectra 
entering ©-(O to the three-dimensional matter power spectrum: 



The matter power spectrum and the distances /kCfX which also 
enter q®(x)> depend on cosmology and can therefore be ex- 
ploited to constrain cosmological parameters. While distances 
and Pss are well known from theory, the probability distri- 
bution of galaxies p^ l) (x) nas to be measured by using ad- 
ditional spectroscopic redshift information ("e.g. iHuterer et alJ 
2006; iMa et alj|200fj; lAbdalla et al J 120071: iBridle & KingT l2007: 
iBernstein & Huterer]|2010h with a certain level of uncertainty. It 
may also be possible to infer s ome additional information from 
the cosmic shear data itself dNewmanl 120081 : ISchneider et al.l 
2006; Zha ng et al.l2 010). The same holds for the slopes of the lu- 
minosity function which can be determined from the survey 
by studying the flux of galaxies close to the magnitude limit. The 
least known quantities in the equations above are the bias terms 
j&g, b\, r g , nj, for which we will thus introduce a very general 
parametrisation in Sect. 13. 31 
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3. Modelling 

In this section we detail the modelling of the terms entering dT7b 
to d26l ). We specify how we parametrise the uncertainty in the 
galaxy redshift distributions, the slope of the luminosity func- 
tion, and the bias terms. Moreover we describe our Fisher matrix 
approach and the way we infer the resulting errors on cosmolog- 
ical parameters. 

3.1. Matter power spectrum & survey characteristics 

As the basis for our analysis we compute matter power spec- 
tra for a spatially flat CDM universe with fiducial parameters 
fi m = 0.25, Q DE = 0.75, and H = 100/iiookm/s/Mpc with 
/iioo = 0.7. We incorporate a variable dark energy model by 
parametrisi ng its equation of state, relati ng pressure pp v. to den- 
sity pde, as dChevallier & Polarsldl200HlLindei1l2003l) 



PDe(z) = \WQ + Waj^-^j PDE(Z) 



(27) 



where the ACDM Universe is chosen as the fiducial model, i.e. 
wo = - 1 and w a = 0. The dar k energy densi ty parameter is then 
given by integrating eq. (3) of iLinderj d2003l) . 

Qde(z) = C^de exp |3 ((w + w a + l)ln(l + z) - W aj^-J} -(28) 
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The primordial power spectrum of matter density fluctuations is 
assumed to be a power law with fiducial slope n s — 1 . We em- 
ploy the fiducial normal isation erg = 0.8. The transfer function of 
Eisenstei n & Hul {1998) is used without baryonic wiggles, com- 
puting the shape parameter with a fiducial value of Qb = 0.05. 

The non-linear corrections to the p ower spectrum are com- 
puted by means of the fit formula by ISrnith et alJ (120031) . We 
account for the influence of dark energy on structure growth by 
modifying the halo mode l fitting routine of I Smith et al.l (2003) 
following the approach of iRefregier et alJ (120081) . We provide a 
summary of this modification in AppendixlDl 

The survey characteristi cs follow the rough specifications 
of a Stage IV experiment ( Alb recht et al.l l2QQ6h such as the 
ESA Euclid satellite mission. To compute the noise proper- 
ties, we assume the maximum extragalactic sky coverage of 
A = 20000 deg 2 and a total number density of galaxies n — 
35 arcmkT 2 . Shape noise is characterised by a total dispersion 
of intrinsic ellipticities of <r e = 0.35. We refer to this survey as 
Euclid-like in th e remainder of this paper. 

According to Smai l et al.l d 19941) we assume an overall num- 
ber of galaxies per unit redshift, per square arcminute 



(29) 



with the galaxy surface density So and B = 1.5. The probabil- 
ity distribution over all the galaxies /?tot(z) is proportional to the 
number density n tot {z). We set z = 0.64, which produces a distri- 
bution with median redshift z mec j = 0.9. The distribution is cut at 
Zmax = 3 and then normalised to unity. For the tomography we 
define photometric bins by dividing the distribution ( 1291 such 
that every bin contains the same number of galaxies. This choice 
is merely for computational convenience and to allow for an easy 
comparison between results with a different number of bins. As 
default we will use N z bm = 10 bins. 

To account for photometric redshift errors, we assume that 
the fraction of catastrophic failures in the assignment of photo- 
metric redshifts is negligible, but include the spread of the true 
redshifts in the bin-wise distributions by writing the conditional 
probability of obtaining a photometric redshift z p h given the true 
redshift z as 

p( -" u, " xp {-i^7TT?j' <30) 

where <x p h denotes the photometric redshift dispersion. The red- 
shift distribution of an individual photo-z bin p^'\z) is then ob- 
tained by integrating (f30b over the bin width and by weight- 
ing the result by the overa ll redshift distribution ( |29l , see 
iJoachimi & Schneider! (120091) for details. We use <x p h = 0.05 as 
our de fault value. This also follows fiducial Model 1 of M a et all 
(2006). 

Since the underlying redshift distributions p®(z) are deter- 
mined by measurement, they are not perfectly known, but intro- 
duce further uncertainty into the analysis. A detailed analysis of 
the dependence of the joint analysis of galaxy shape and num- 
ber density information on redshift parameters, and also the po- 
tential of calibr ating these errors i n ternally, will be invest igated 
elsewhere (e.g. IZhang et al.ll2010l) . iBridle & Kingl d2007l) have 
undertaken a more detailed study of the effect of redshift errors 
in the case of ellipticity correlations only. For the purpose of this 
work we assume that the value of <r p h is unknown, i.e. we use it 
as a single, global parameter to account for the uncertainty in the 
redshift distributions. We employ a wide Gaussian prior on cr ph 
of 10 for reasons of numerical stability. 



3.2. Galaxy luminosity function 

In order to calculate power spectra which include the lensing 
magnification signal, we need to model the slope of the cumu- 
lative galaxy luminosity function at the magnitude limit of the 
galaxy number density catalogue. In Appendix ICl we extend ob- 
servational results for the normalisatio n and redshift s caling of 
the galaxy redshift distribution d29]> bv lBlake & Bridld d2005ri to 
provide a fitting formula for the luminosity function slope as a 
function of redshift and survey magnitude limit. 

We use the fit given by dC.2l > and dC.3l ) with the parameters 
listed in Table IC.ll to compute the slope of the luminosity func- 
tion at ri; m = 24. The discussion in this work applies to ground- 
based surveys because the COMBO- 17 luminosity functions are 
calculated for the SDSS r filter as observed from the ground. A 
space mission to a depth of r\i m - 24 will have a different lu- 
minosity function slope, corresponding more closely to a deeper 
ground based survey, depending on the resolution of the space- 
based survey. We use results for r\i m = 24 throughout this paper 
for both ground and space surveys. We note that from Fig. IC.ll 
top panel, the slope of the luminosity function is changed little 
on increasing the survey depth beyond n; m = 24. 

The fiducial slope in a photo-z bin i is defined as or® = 



cz(z , , ni m = 24), where z w , is the median redshift of bin i, 



_(0 



med 



see Appendix lAl We assume ar' is also measured from the sur- 
vey itself, and therefore adds another source of uncertainty to the 
analysis which we account for by setting a w to be free parame- 
ters for all i = 1 , .. , A^ z bin- Again we apply a wide Gaussian prior 
of 10 on every slope parameter, which does not have a mea- 
surable influence on our results, but merely ensures numerical 
stability. 

3.3. Galaxy and intrinsic alignment bias 

As already outlined in Sect.|2] the bias terms bx = \b g , bj, r g , 
encoding the galaxy bias and intrinsic alignments are the least 
accurately known contributions to dTD l to d26l >. We parametrise 
each o f these terms on a grid in k and z, following lBridle & Kind 
(2007) whose ansatz is in tur n similar to the recom men dations by 
the Da rk Energy Task Force (Albr echt et al.f 2006) and Bernstein 
(2009). Every bias term is assumed to vary around a fiducial 
functional form fr^ ase as 



bx(Kx) = A x Qx(k,z(x)) 



(31) 



where Qx(k,z(x)) is an unknown two-dimensional function 
which comprises the aforementioned grid, and where Ax denotes 
an additional free overall amplitude. We use Nk bins in k and Nz 
bins in redshift for each bias term and linearly interpolate in the 
logarithms of Qx , k, and 1 + z, so that Qx is given by 

In Qx (k, z) = K(k) Zj(z) Bf } + [l - ZjQO B? i+l)j (32) 

+*,(*) [l-Zj<d] Bf u+1) + [l - K(k)] [l-Zj®] Bf i+l)(j+1) 
for kj < k < k i+ i and Zj <Z< Zj+i, where we defined 
ln(fe) - Info) 



K,(k) = 
Zj(z) = 



lnfo +1 )-lnfo) 

ln(l +z)-ln(l + Zj ) 
ln(l + zj+i) - ln(l +zj) 



(33) 



The free parameters are the grid nodes Bfj. Since Qx is a multi- 
plicative function, one reproduces the base model Z?^ ase if Ax = 1 
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Table 2. Overview on the total number of nuisance parameters 
used for different setups. 



i.e. they are fully determined by the base models. We set 



b h ™{k, X ) = 1 



N K xN z 


66 


nn 


all 


2x2 


11 


21 


31 


2x4 


19 


29 


47 


10 x 2 


43 


53 


95 


10x4 


83 


93 


175 


3x3 


21 


31 


51 


5x5 


53 


63 


115 


7x7 


101 


111 


211 



Notes. The label ee corresponds to using ellipticity correlations only 
as the observables. Likewise, nn corresponds to using galaxy number 
density correlations only, and 'all ' to using all available correlations. 
For all entries we have assumed ten photometric redshift bins used for 
the tomography, A^in = 10. 



and if all Bfj vanish. The effect of this parametrisation on the ob- 
servable projected power spectra is illustrated in lBridle & Kind 
(2007) for the case of intrinsic alignments. 

The indices in (1321 run from i = 0, . . , Nk + 1 and 
I . .. ,Nz + 1. We fix all parameters at the edge of the grid by 
setting the parameters with indices i,j = 0; i = Nk + 1 or 
j — Nz + 1 to Bfj = 0, so that we have Nk X Nz free grid param- 
eters per bias term. We place the lowest and highest grid nodes 
at the limits of our integration ranges, so ko = 3.3 • 10 _7 ftMpc _1 
and kN K+ \ = 3.3 • 10 4 ftMpc _1 in k, and in the redshift dimen- 
sion zo = and zn z +\ = 19. The grid nodes, which are free to 
vary, are log-linearly spaced in a smaller range, respectively. We 
use k\ - 10~ 3 /i Mpc -1 and k^ K = 2ft Mpc, and for the redshift 
range z\ = Zo and zn z = 3. In the special case of Nk — 1 we 
position the only free parameter in the k dimension at the centre 
between k\ and k^ K , and proceed likewise for redshifts. 

It is important to note that while Bri dle & Kind (l2007h limit 
the flexible grid parametrisation to the non-linear regime of the 
power spectra, we attempt to cover all k ranges which substan- 
tially contribute to the observable power spectra. As we fix the 
grid values on the edges, the overall scaling of the bias terms 
is not free, so that we use the amplitude Ax as a further vary- 
ing parameter throughout. To all bias term parameters we add 
a very wide Gaussian prior of standard deviation 50 to ensure 
numerical stability. Together with the global uncertainty on the 
redshift distributions, expressed in terms of <x ph , and the val- 
ues of the slope of the galaxy luminosity function per photo- 
z bin, we obtain a large number of nuisance parameters that 
we determine simultaneously with the cosmological parameters 
of interest. For later reference, we have summarised the total 
number of nuisance parameters for different setups in Table [2] 
While our parametrisation is fairly general and should capture 
most of the variability, it is of course possible that the bias 
terms depend on more parameters than k and z. For instance, 
it is well known that both intrinsic alignments and galaxy bias 
are a function of galaxy colour and luminosity which could 
be incorporated into our approach in the future. For observa- 
tional constraints on this e ffect using intrinsic alignments see 
iMandelbaum et al.l (120061): iHirata et'alJ (120071) and in galax y 



biasing see McCracken et al. (2008); Swanson et al. (2008); 



Simo n et al.l (120091) : ICresswell & Percivall (120091) : IWang et all 
d2007l) for recent examples. 

To compute the fiducial models for all the power spectra en- 
tering the observables (HJi to (O, we set Ax — 1 and all Bf. - 0, 



g 

base 



(34) 



(k,x) 



P(z) 



D(z)(l+z) 



r b r(k,x) = 1 • 

With the choice for b^ ase and r^ se we reproduce the non-linear 
modifi cation of the linear alignment model by iBridle & Kind 
(12007b . Lacking solid physical motivati on, it is yet in agreement 
wit h current observation al evidence (Mandelba um et al.l 12006; 
see IBridle & Kindl2007l for a c ompa rison) and the halo model 
studies by ISchneider & Bridiel (l2010l). It is based on the linear 
alignment model (Hirata & S el iakl 120041) which is expected to 
provide a good description of intrinsic alignments on the largest 
scales. We assume the galaxy bias to be of order unity for our 
fiducial model, and set r h ™ e (k,x) — 1- Note that most investi- 
gations of galaxy clustering consider much less flexibility in the 
galaxy bias. 

In Fig.[TJ we plot the fiducial angular power spectra of all 
considered signals for different combinations of photo-z bins. 
Note that the ellipticity-number density cross-correlations are 
not symmetric under exchanging the photometric redshift bins. 
Hence, in this figure we treat ne and en correlations, as well as 
all signals contributing to them, separately, keeping i < j for all 

The GG signal shows the usual behaviour of moderate in- 
crease with increasing redshift of the contributing photo-z-bins. 
The redshift scaling of the IG term is similar, but peaks when the 
source galaxies are at high redshift in the background (large j), 
while the galaxies that are intrinsically aligned are at low red- 
shift (small i). For the model used here the IG contribution can 
even surpass the shear signal in this case. Due to the narrow ker- 
nel containing p^'Xx) p (x)> see ©> the II signal is strong in 
the auto-correlations i = j, but drops off quickly as soon as the 
overlap of the redshift distributions decreases. 

Due to the similar kernel, the scaling of the galaxy cluster- 
ing contribution (gg) resembles the II term, but gg constitutes 
a much stronger signal. Lensing magnification (mm) adds the 
largest fraction of the galaxy number correlations at the high- 
est redshifts, showing a slightly stronger redshift scaling than 
GG. However, the mm term always remains subdominant with 
respect to signals with a contribution from galaxy clustering; 
even for widely separated galaxy redshift distributions, say i - 1 
and j - 9, and the gm cross-term is considerably stronger than 
mm. Such contributions might be a serious obstacle for prob- 
ing cosm ology with the l e nsing magnification s ignal as pro - 
posed by Broadhur st et all (1 1995b : IZhang &Penl (120051 [2006); 
van Waerbeke ( 2010). Yet in our approach, where the galaxy bias 
is taken into account and parametrised, the magnification signal 
yields a valuable contribution to the galaxy number correlations, 
which helps constraining the cosmological model. 

The signals within the number density-ellipticity cross- 
correlations are not symmetric when swapping the photo-z bins. 
When the contribution by number density fluctuations stems 
from the foreground, the gG signal is strong, in particular if the 
photo-z bins are far apart in redshift, whereas the 'Gg' (in the 
notation of Fig.[TJi drops off fast if i < j because the shear sig- 
nal of foreground galaxies is not correlated with the clustering 
of galaxies at much higher redshift. The mG, Gm, and GG sig- 
nals differ only by the term including the slope of the luminosity 
function and thus have similar amplitudes. Correlations between 
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Fig. 1. Fiducial power spectra for all considered correlations. The upper right panels depict the contributions to ee (in black) and 
nn (in magenta) correlations. The lower left panels show the contributions to correlations between number density fluctuations and 
ellipticity. Since we only show correlations C^{€) with i < j, we make in this plot a distinction between ne (in red; number density 
contribution in the foreground, e.g. gG) and en (in blue; number density contribution in the background, e.g. Gg) correlations. In 
each sub-panel a different tomographic redshift bin correlation is shown. For clarity only odd bins are displayed. In the upper right 
panels the usual cosmic shear signal (GG) is shown as a black solid lines; the intrinsic alignment GI term is shown by the black 
dashed lines; the intrinsic alignment II term is shown by the dotted black line; the usual galaxy clustering signal (gg) is shown by 
the magenta solid line; the cross correlation between galaxy clustering and lensing magnification (gm) is shown by the magenta 
dashed line; the lensing magnification correlation functions (mm) are shown by the magenta dotted line. In the lower left panels the 
solid blue line shows the correlation between lensing shear and galaxy clustering (Gg); the blue dashed line shows the correlation 
between lensing shear and lensing magnification (gm); the blue dot-dashed line shows the correlation between intrinsic alignment 
and galaxy clustering (Ig or equivalently gl); the red solid line shows the correlation between galaxy clustering and lensing shear 
(gG), which is equivalent to the blue solid line with redshift bin indices i and j reversed; similarly the red dashed line shows the 
correlation between lensing magnification and lensing shear (mG), for cases where the magnification occurs at lower redshift than 
the shear (i < j); finally the dotted line shows the correlation between lensing magnification and intrinsic alignment (ml). 
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intrinsic alignment and magnification (ml and Im) are subdom- 
inant throughout, obtaining their largest amplitudes if intrinsic 
alignments at low redshifts are combined with the magnifica- 
tion signal from galaxies far in the background, i.e. in the figure 
for Im at i = 1 and j = 9. Finally, the symmetric Ig term is 
the largest contribution for auto-correlations of number density- 
ellipticity observables, but decreases quickly in the cross terms, 
again due to the kernel p®(x) p^(x) since we have assumed the 
photometric redshift errors are reasonably well behaved, with- 
out catastrophic outliers. Note that we have plotted the absolute 
values of the power spectra in Fig.Q]and that the correlations IG, 
gm, Gm, mG, and gl are negative. 

3.4. Parameter constraints 

We determine constraints on our parameters using a Fisher ma- 
trix analysis. To account for the errors and correlations of our 
observables, we compute covariances of the power spectra © 
to © in the Gaussi an approxi mation, extending the results of 
Joachimi et al. (2008); see also Hu & Jai 3 (|2004). If we denote 
the difference between estimator and its ensemble average by 
iSC^{€), one can write for the covariance 



(ACg>(0 AC^V')) = Ste 



2n 

AIM 



(35) 



where A is the survey size and A£ the width of the corresponding 
angular frequency bin. As the Rronecker symbol dec indicates, 
the covariance is diagonal in I in the Gaussian limit, which keeps 
the computation and inversion of ( T35l l tractable. The subscripts 
[<x,p, y, 5} can be either e or n, where C^(f) = C { ^{€) holds. To 
account for the shot and shape noise contributions induced by 
the random terms in ([TJ and ©, we have defined 



the second term given by 

J2 



N' 



1 



(37) 



m 

N™ = 5 U — fora=/J = « 
=0 for a # p . 

Here cr\ denotes the total intrinsic ellipticity dispersion, and n w 
is the average galaxy number density per steradian in photo-z 
bin i. 

Combining the observable power spectra, we compose the 
total data vector 

D{€) = jc^OD, ..,C { ^ N ^\i),Cg l \i), .., (38) 

for every angular frequency considered. The corresponding co- 
variance, again for every I, reads 



Cov(0 



Cov^SV) 


Cov^V) 


Cov^W 


Cov^(0 


Cov» 


CovJ«S?(0 









(39) 



with the block matrices given by d35l >. The number of galaxy 
ellipticity (ee) and number density (nri) observables entering 
D({) is A'zbin (A^zbin + l)/2, respectively, while there are N^ bi 
ellipticity-number density cross terms (ne), which are not sym- 
metric. In the analysis that follows we will also consider ee and 
nn correlations only. In these cases the covariance ( |39l is re- 
duced accordingly to its upper left or lower right block. 

For reasons of computational time the total number of pa- 
rameters that we can consider is limited to a few hundred. As 
a consequence the k dependence of the galaxy bias can not be 
parametrised by more than about ten parameters per redshift 
grid node. This number might not provide enough freedom in 
b g (and r g ) to represent a sufficiently general set of functional 
forms, which inadvertently may cause strong constraints on cos- 
mological parameters due to the strong signal of galaxy cluster- 
ing. Hence, we follow existing studies of galaxy clustering by 
discarding the clustering contribution in the non-linear regime 
where the signal is largest and the form of the galaxy bias most 
un certain. 

iRassat et al.l d2008l) calculated wave vectors fc™ x as a func- 
tion of redshift at which the three-dimensional power spectrum 
has to be cut off to avoid number density fluctuations above a cer- 
tain threshold, used as an indicator for non-linearity. Since we do 
still have a fairly general parametrisation of the bias terms asso- 
ciated with galaxy bias, we can afford to include the mildl y non- 
linear regime into our analysis. Consulting Fig. 2 of Rassa t et al.l 
d2008h . we choose a simple linear parametrisation of the form 



kZ x (z)~ 0.132 zh Mpc- 



(40) 



Min 

This relation roughly coincides with the fiducial curve in the fig- 
ure, producing slightly more conservative cuts at low redshifts. 

We do not cut the three-dimensional power spectrum in 
fc-space, but instead exclude projected power spectra above a 
threshold angular frequency from the Fisher matrix analysis. 
This maximum angular frequency is computed via 

(36) e=0&AM. (41) 

where we choose as a characteristic redshift of bin i the median 
redshift z w .. Hence, we obtain a cut-off £ for every photo-z bin. 



We choose that ee correlations are not at all affected by this cut- 
off because they are not dominated by terms involving galaxy 
bias. We impose {^i on ne correlations, where i is the photo-z 
bin from which the number density signal stems. For observables 
C„n{€) we use the cut-off calculated for bin j. Note that, due to 
the fast drop-off of the galaxy clustering signal with increasingly 
different median redshifts of bins i and j, the more optimistic 
choice of bin j over i in the latter case should not influence our 
results decisively. 

To compute the Fisher matrix, we use Ne = 50 angular 
frequency bins, spaced logarithmically between l m \ n = 10 and 
Anax = 3000, the latter value being a conservative maximum 
for future surveys. We assume that the covariance of the power 
spectra is independent of the cosmological parameters, so that 
it does not contribute to the c onstraints. Then the Fisher matrix 
reads (e.g. Tegmark et al l 19971) 



Nd W f m «(m,n) 

Z Z 



Cov m J£) 



(42) 



dp, mn " ' d Pv 

where Nd is the dimension of T>(€). The first summation in d42l) 
runs over all A^™ ax usable angular frequency bins. The number 
of usable bins depends on the type of correlation and is deter- 
mined by the cut-off angular frequency as described above. For 



9 



B. Joachimi and S. L. Bridle: Simultaneous measurement of cosmic shear and galaxy number density correlations 



Table 3. Overview on the cuts in angular frequency. 



the error ellipsoid in parameter space, which is measured by the 
determinant of the Fisher matrix. Hence, we define 



z-bin 




c& 

'max 


Annas 


t 
i 


O.ij 


A c 
43 


13 


2 


0.51 


103 


20 


3 


0.63 


145 


23 


4 


0.74 


206 


26 


5 


0.85 


260 


28 


6 


0.96 


329 


30 


7 


1.08 


369 


31 


8 


1.23 


466 


33 


9 


1.42 


588 


35 


10 


1.79 


937 


39 



Notes. Results are shown for the default set of parameters with A'zbin = 
10 and iTph = 0.05, using J40t and J4 It . For each photo-z bin the median 
redshift z mc(i , the maximum angular frequency £f mx , and the correspond- 
ing number of usable angular frequency bins N™ ax (out of the total of 
50) is given. 



-<<*> 



/>? <<") 



and the number 



the default setup we have summarised z med - 
of usable bins for every photo-z bin in Table [3] The derivatives 
in (l42l are taken with respect to the elements of the parameter 
vector 

p = |Q m ,cr 8 ,/j 100 ,n s ,a b , w Q , w a ; (43) 
tr p h,u; ,..,u ,D n , .. ,o NkNz 

The first seven entries of p correspond to the cosmological pa- 
rameters that we are interested in, while the remaining nuisance 
parameters account for the uncertainty in the galaxy redshift dis- 
tribution, the slope of the galaxy luminosity function, the intrin- 
sic alignments, and the galaxy bias, see Table [2] In summary, 
we use a maximum of Ne{2N^ b - m + N z u D ) observables (actually 
significantly less due to the /'-cuts of the galaxy number density 
signals) to measure a total of 4-NzNK+N z u n + 8 parameters. Note 
that since we have referred all signals contributing to the observ- 
ables © to (0 to the matter power spectrum, they all constrain 
the set of cosmological parameters; none of them is fixed when 
calculating d42l >. 

The minimum variance bound of the error on a parameter 
Pfj, if determined simultaneously with all other parameters, is 

given by crip^) = y/iF' 1 )^. This error provides us with a lower 
bound on the marginalised 1 <x-error on p^. To assess the statisti- 
cal power of the survey by means of a single number, we use the 
figure of merit ( FoM) suggested by th e Dark Energy Task Force 
(DETF) Report dAlbrecht et al.ll2006l) . 



FoM 



1 



DETF 



(44) 



where the subscript 'wow a ' denotes the 2x2 sub-matrix of the 
inverse Fisher matrix that corresponds to the entries belonging 
to the two dark energy parameters. Note that different prefactors 
for ( l44t a re used in the lit erature. To allow for direct comparison 
with Bridl e & Kingl (2007), divide our findings for the FoMdetf 
by four. 

While (PPfl l is restricted to the quality of constraints on dark 
energy, we also seek to consider the errors on all cosmological 
parameters of interest. We are interested in the total volume of 



FoMtot = In 



1 



det(F-!) c 



(45) 



where only the sub-matrix of the inverse Fisher matrix that cor- 
responds to the seven cosmological parameters under investiga- 
tion is used in the determinant, as indicated by the subscript. The 
determinant of the inverse is computed in order to take the effect 
of marginalising over nuisance parameters into account. 

4. Results 

Based on the Fisher matrix formalism described in the forego- 
ing section, we will now analyse the performance of a cosmo- 
logical galaxy survey with combined number density and shear 
information. We are going to investigate the residual informa- 
tion content in the data after marginalising over models of the 
intrinsic alignments and the galaxy bias with varying degrees 
of freedom. Furthermore, we will study the dependence of our 
FoM on the number of photo-z bins and the width of the bin- 
wise redshift distributions as well as on the priors imposed on the 
nuisance parameters. The information contained in the individ- 
ual signals and their susceptibility to the nuisance parameters is 
also assessed. Throughout this section we use the default survey 
characteristics and parameter values unless specified otherwise. 

4.1. Dependence on intrinsic alignments and galaxy bias 

The four bias terms \b g , b\, r % , rA each comprise Nk X Nz + 1 
nuisance parameters. The galaxy ellipticity (ee) power spectra 
only contain the two intrinsic alignment bias terms, whereas the 
number density correlations (nn) are only affected by the galaxy 
bias terms. The cross-correlations between ellipticity and num- 
ber density link those signals which depend on both galaxy bias 
and intrinsic alignments and thus allo w for their inte rnal cross- 
calibration. An example is the study bv lZhangl (l2010) which in- 
vestigates the interrelations between the IG, gl, and gg terms. 

In Fig. [2] both figures of merit are plotted as a function of 
Nk and Nz, respectively. If we restrict our analysis to ee corre- 
lations on ly, our setup is s imilar to the most realistic setup con- 
sidered in Bridle & King] (120071) . We reproduce their result that 
the FoMdetf drops as a function of Nk, dropping most sharply 
at small Nk and then levelling off. It falls significantly below the 
reference value of the FoMdetf computed for a pure GG signal. 
However, while iBridle & King d2007l) find a decrease by about a 
factor of 2 for Nk =10 and Nz = 2 compared to the lensing only 
case, our FoMdetf decreases by more than a fa ctor of 4. This 
discrepancy can be traced back to the fact that Bridle & King 
J2007h limit the nuisance parameter gridding to the nonlinear 
regime in £-space. We additionally plot the figures of merit as a 
function of Nz, finding that the ee results do in fact flatten as Nz 
reaches high values around 10. We find a very similar behaviour 
in terms of the FoMtot- The dark energy parameters suffer more 
than other parameters from the uncertainty of the bias terms in 
the redshift direction, as the FoMdetf deteriorates faster than the 
FoMtot as a function of Nz- 

With our default settings the pure nn correlations con- 
strain cosmology much more weakly than ellipticity correla- 
tions. Recall that our galaxy clustering analysis uses a much 
more flexible bias parametrisation than most other work. Using a 
small number of nuisance parameters we get a FoMdetf which 
is marginally above unity. This result is of the same order of 
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Fig. 2. Left panels: Figures of merit as a function of the number of free parameters as a function of wave vector Nk in the bias 
terms For each line type, the upper curve is obtained for a number of free bias parameters as a function of redshift Nz = 2, the 
lower is for Nz = 4. Right panels: Same as on the left, but as a function of Nz, i.e. the number of redshift parameters in the bias 
terms. The upper curves for each set correspond now to Nk = 2 and the lower ones to Nk = 4, respectively. Upper panels: Figure of 
merit taking into account the full cosmological parameter space, FoMtot, see ( 1431 ). Lower panels: Dark Energy figure of merit from 
the Dark Energy Task Force FoMdetf, see d44b . Dashed curves correspond to results using galaxy ellipticity correlations (ee) only, 
dotted black curves to galaxy number density correlations (nn) only, and solid black curves to results using all correlations (ee, nn 
and en). The grey dotted lines show results for nn correlations without imposing cuts in angular frequency. The constant grey line 
marks the FoM computed for the pure lensing, i.e. GG, signal, assuming intrinsic alignments do not exist. In addition we show the 
resulting figures of merit when using our most flexible parametrisation with Nk — Nz — 7 as filled symbols. Circles correspond to 
ee, triangles to nn, and diamonds to all correlations. 



magnitude as the findings of iRassat et al.l d2008l) who deter- 
mine FoMdetf = 4.2 for a spectroscopic, space-ba sed survey in 
the sp herical harmonics approach. Contrary to us, Rassa t et alJ 
(2008) use a less flexible bias parametrisation and include BAO 
features, but neglect magnification effects in their study. The de- 
crease in both figures of merit for a larger number of nuisance 
parameters is more pronounced than for ee correlations, in par- 
ticular as a function of Nz- 

We have also shown the results obtained without the addi- 
tional angular frequency cuts in Fig.[2](i.e. the maximum angular 
frequency is ^ max = 3000 for all angular power spectra). Because 
the full galaxy clustering signal is strong and a comparatively di- 
rect probe of the dark matter power spectrum if the galaxy bias is 
well known, the values of the figures of merit are much higher in 
this case and even surpass the lensing-only level if few nuisance 
parameters are used. This corresponds to the case where we have 
a reasonably good understanding of galaxy biasing over a large 
range of scales, into the deeply non-linear regime. For larger Nk 
and Nz the decrease in FoM is considerable, but weaker than 
with ^-cuts. 

In the limit of a large number of nuisance parameters we ex- 
pect that the galaxy bias is parametrised in a sufficiently flexible 
way, such that the curves with and without /"-cuts should coin- 



cide or at least be of the same order of magnitude. Since this is 
not the case, and since there is no obvious sign of a lower bound- 
ary that the figures of merit are approaching, we hypothesise that 
this limit is only achieved for a very large, numerically and com- 
putationally prohibitive, number of nuisance parameters. Thus, 
the use of the additional angular frequency cuts in our default 
analysis is the most practical way to take into account our lack 
of understanding of galaxy biasing on small scales. 

The simultaneous use of all available observables signif- 
icantly boosts the parameter constraints, due to the addition 
of extra information from the survey, and due to breaking de- 
generacies between cosmological and nuisance parameters. The 
FoMdetf is up to a factor of about 50 higher than the lensing- 
only value. Both figures of merit decrease for larger Nk and Nz, 
attaining mostly shallow slopes at Nk,Nz ~ 10. Considering the 
most flexible model in each panel of Fig.|2l the FoMtot remains 
above, but close to the value for a pure GG signal while the 
FoMdetf degrades slightly below the lensing-only value. In con- 
trast to the setup with ee or nn correlations only, the decrease is 
less pronounced with Nz than Nk here, i.e. the cross-calibration 
between all correlations can partially compensate the loss in red- 
shift information that affected the dark energy parameters in the 
former cases. 
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Table 4. Marginalised parameter errors and the two Figures of Merit for a survey with ten photometric redshift bins used for 
tomography N z b; n = 10 and a photometric redshift uncertainty parameter <x p h = 0.05. 



par. 




Euclid-like survey 






DES-like survey 






ee 


nn 


all 


lensing 


ee 


nn 


all 


lensing 


i l in 


0.0228 


0.1530 


0.0079 


0.0044 


0.0636 


0.3620 


0.0233 


0.0120 


erg 


0.0295 


0.1800 


0.0110 


0.0065 


0.0827 


0.5941 


0.0335 


0.0180 


hioo 


0.2072 


0.4273 


0.0647 


0.1099 


0.5447 


1.1092 


0.1918 


0.2837 




0.0721 


0.2218 


0.0312 


0.0250 


0.2117 


1.0130 


0.0973 


0.0661 


fib 


0.0224 


0.0329 


0.0063 


0.0133 


0.0557 


0.0771 


0.0166 


0.0341 


>>'() 


0.3242 


1.6404 


0.0939 


0.0579 


0.9064 


4.3967 


0.2895 


0.1631 


w a 


1.1995 


5.6187 


0.3051 


0.2014 


3.4754 


18.2522 


0.9865 


0.5728 


FoMtot 


48.30 


24.80 


67.13 


66.12 


33.64 


7.70 


51.87 


52.29 


FoM DETF 


9.12 


0.28 


102.92 


182.96 


1.08 


0.04 


10.96 


24.92 



Notes. We have used the configuration with the maximum set of nuisance parameters that we consider, i.e. N K = N z = 7. The cuts in angular 
frequency for number density correlations have been applied to restrict the analysis to the linear regime. Shown are the results using shear- 
shear correlations ee only, observed galaxy clustering correlations nn only, all correlations, and the lensing-only signal (assuming no intrinsic 
alignments). The left-hand part of the table uses the standard Euclid-like survey parameters, i.e. <x f = 0.35, n = 35 arcrrmT 2 , and the survey area 
A = 20000 deg 2 . The right-hand part of the table uses parameters for a DES-like survey, i.e. <x f = 0.23, n = 10arcmin~ 2 , and the survey area 
A = 5000 deg 2 . 



The most flexible configuration we consider has Nk — Nz = 
7 which corresponds to 200 nuisance parameters within the four 
bias terms. This limit is not inherent to our method, but is merely 
set for computational practicality. Both the intrinsic alignments 
and the galaxy bias have a physical origin and therefore are ex- 
pected to produce smooth signals, which should not oscillate 
strongly or feature sharp peaks. Hence on the default angular 
scales used this model with 100 free parameters for each of in- 
trinsic alignments and galaxy bias should yield a fairly general 
representation of the signals if one can rely on coarse prior in- 
formation on the fiducial base model. 

The figures of merit and the individual parameter errors for 
the most flexible model are given in Table [U for the pure GG 
signal, ee, nn, and all correlations. Compared to lensing alone, 
the FoMdetf decreases by about a factor of 20 for ee correla- 
tions, which means that cosmic shear is severely affected if one 
assumes very little prior knowledge about intrinsic alignments. 
There are hardly any dark energy constraints for nn correlations 
only, given this freedom in the galaxy bias. However, using all 
available correlation simultaneously, we can recover just over 
half of the pure cosmic shear DETF figure of merit. We find that 
the total error volume of the cosmological parameters, i.e. the 
FoMtot, is the same using all the correlation information as for 
the lensing-only case. Looking at the marginalised parameter er- 
rors, wo and w a suffer particularly strongly in the ee only and nn 
only case. Less information about the dark energy parameters is 
lost if all correlations are used, and the constraints on /iioo and 
Qb improve over the lensing-only case. In Fig. [3] we show in ad- 
dition the marginalised 1 cr-error ellipses of all possible pairs of 
cosmological parameters for the same setup. The degeneracy di- 
rections are fairly similar for the usual lensing only case and the 
results with all correlations presented in this paper. 

We repeat this analysis for a nearer term Stage-III-like sur- 
vey such as the Dark Energy Survey (DES). To this end, we use 
a survey size of A — 5000 deg 2 , a median redshift of the over- 
all redshift distribution of z me d = 0.8, a total number density of 
galaxies of n — 10 arcmin~ 2 , and a total ellipticity dispersion of 
<r e = 0.23, keeping all other parameters at their default values. 
We refer to this survey as DES-like for the rest of the paper. 

Note that different groups use different conventions for quot- 
ing the number density of galaxies and the ellipticity dispersion. 
The Euclid team define <x e = 0.35, which is motivated by con- 



sidering it to be the typical uncertainty on the shear taking into 
account both shape noise and measurement errors. The effective 
number density has then been estimated rel ative to this value an d 
found to be between 30 and 40 arcmin ~ 2 dLaureiis et alJ l2009). 
Because of the choice of cr e , this number density value is ap- 
proximately equal to the actual number of usable galaxies in the 
image. The DES team start by defining <r e - 0.23, motivated by 
considering the shape noise alone. T he effective numbe r density 
is then found to be n — 10 arcmin 2 dAnnis et alj|2005l) . Due to 
this choice of cr e , the number density value is smaller than the 
number density of usable galaxies in the image, and represents 
the effective number density of galaxies with negligible mea- 
surement errors. The different definitions of cr e also explain the 
low number density quoted for DES in comparison with shal- 
lower surveys such as the CHFTLS. 

In this work we choose to take the number density and el- 
lipticity dispersion values as provided by the collaborations and 
do not adjust them, which would be beyond the scope of this pa- 
per. In any case the physically meaningful quantity is the ratio 
cr e / V» which determines the uncertainties on the power spectra, 
see ( l37l l. Thus, as long as we use consistent definitions of n and 
cr e , we obtain the correct shape noise contributions to the cosmic 
shear covariance. 

Our findings are also shown in Table |4] For the DES-like 
survey the FoMtot for all correlations is again about the same 
as in the lensing-only case, and the ratios of FoMdetf values are 
slightly smaller than for the Stage-IV-like survey. The FoMdetf 
values for ee, nn, and all possible correlations are about a factor 
of ten larger from the Stage-IV-like cosmic shear survey than 
the Stage-III-like survey , offering significan t benefit beyond the 
minimal requirement of Alb recht et al.l (120061) . 

4.2. Dependence on characteristics of the redshift 
distribution 

It is well known that dividing the galaxy sample into several 
redshift photo-z bins greatly improves constraints from cos- 
mic shear, but due to the broad lensing kernel (TTOb there is lit- 
tle benefi t in having more d i stributions than three to five (e.g. 
lHulll999t ISimon et all 12004 iMa etai1l2006h . This result does 
not hold true anymore if one aims at controlling the intrin- 
sic alignment contamination in the cosmic shear signal, mainly 
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Fig. 3. 1 cr-contours for all pairs of cosmological parameters considered, marginalised over all other parameters. We have used a 
photometric redshift uncertainty parameter value <x ph = 0.05, ten photometric redshift bins for tomography A^ z bi n = 10, and the 
most flexible intrinsic alignment and bias model considered in this paper, with over two hundred free parameters (Ng = 7, and 
Nz = 7). Orange (light hatched) confidence regions result from using galaxy number density correlations (nri) (excluding the non- 
linear regime) only, red (dark hatched) regions use ellipticity correlations (ee) alone, and blue (filled) regions correspond to using 
all available information including density-ellipticity cross-correlations. For reference, the contours obtained from a pure lensing 
signal are shown as black lines. Flat priors on cosmological parameters have been applied. 



manifest via its characteristic redshift dependence. Using both 
marginali sation ([Bridle & Kind 120 07) and parameter-free ap- 
proaches dJoachimi & Schneiderl20 08). one finds that the figures 
of merit only start to stabilise when using ten redshift distribu- 
tions or more. 

In Fig.|4]our figures of merit are plotted as a function of the 
number of photo-z bins, using Nk — Nz — 5 nuisance parameters 
per bias term. In agreement with the findings mentioned above 
both the FoMtot and the FoMdetf in the case of a pure cosmic 
shear signal become approximately constant for N z \>m ^ 5, and 



increasing A^ z bi n beyond three has little effect. For ee correlations 
the gain in FoM is considerable up to A^ z bi n ~ 7; for larger N z hm 
the curves rise only slowly. 

Since nn correlations are more localised because of the com- 
pact kernel of the dominating galaxy clustering (gg) signal, in- 
creasing A^ z bin proves beneficial for these observables up to the 
maximum number of photo-z bins we have considered, although 
the bins feature an increasing overlap of their corresponding red- 
shift distributions as we keep <x p h = 0.05 fixed. As one would 
expect, we obtain an intermediate scaling with A^ z bi n for the 
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Fig. 4. Upper panel: The Figure of Merit for all cosmologi- 
cal parameter space FoMtot as a function of the number of 
photo-z bins used for tomography iV z bi n , shown for ee (dashed 
line), nn (dotted line), and all (solid) correlations. The grey 
line corresponds to results for lensing only (GG). Throughout, 
Nk = Nz - 5 nuisance parameters for the bias terms are used. 
These results are obtained for the standard set of parameters and 
cr ph = 0.05. Lower panel: Same as above, but in terms of the 
dark energy figure of merit FoMqetf- 

complete set of available correlations. Our fiducial choice of 
A^zbin = 10 is beyond the regime of strongly varying figures of 
merit at small iV z bi n , but the further increase in FoM is more pro- 
nounced than for a pure lensing signal or ee correlations only 
with FoMqetf rising by an additional 80 % on increasing the 
number of photometric redshift bins from 10 to 20. 

Figure|5]shows the figures of merit as a function of the photo- 
z dispersion, normalised to the value at <x p h = 0.02. Since the 
FoMtot is a logarithmic quantity, we compute differences rather 
than ratios of the figure of merit, i.e. 

r Fo M = FoM DE TF(o-ph)/FoM DET F(o-ph = 0.02) (46) 
d FoM = FoMxot(°" P h = 0.02) - FoMjo-Ho-ph) (47) 




[det(F-i) cos J (crp h = 0.02) 



Hence the difference t/poM is directly related to the change in 
volume of the error ellipsoid spanned by the set of cosmological 
parameters. We have returned to our default value of A^ z bi n = 10. 

The pure lensing signal needs merely coarse redshift infor- 
mation to attain its full statistical power, and hence its FoM 
hardly suffers from the increasing spread in the redshift distri- 
butions. However, as redshift information is vital to account for 
intrinsic alignments, the figures of merit for the cosmic shear 
signal after marginalisation over intrinsic alignments decrease 
considerably by more than 80 % in terms of the FoMdetf on 



Fig. 5. Upper panel: The difference 6?foM> defined in d47l >. as a 
function of the photo-z dispersion cr ph , shown for ee correlations 
(dashed line), all correlations (black solid line), and the lensing 
only signal (grey solid line). Throughout, nuisance parameters 
Nk = Nz = 5 are used. These results are obtained for the stan- 
dard set of parameters and N z \,i n = 10. Lower panel: Same as 
above, but in terms of the ratio rp M, given in d46l) . 



changing from cr n h = 0.02 to cr p b = 0.1, which is in line with 
iBridle & Kingl (2007). For the same change in cr ph the FoMtot 
is reduced by about 12. It is interesting to note that the depen- 
dence of the FoMtot on cr ph is close to linear for both lensing 
and ee signal. 

Adding galaxy number density information largely alleviates 
the information loss. The degradation in FoM is only slightly 
stronger than for the pure lensing signals in the expected regime 
of high-quality photo-z information with <x ph < 0.06. For larger 
(Tph the increasing overlap of the redshift distributions causes 
both figures of merit to decrease further. We have to add the 
caveat that in this investigation we consider only a single param- 
eter that accounts for the uncertainty in the shape of the redshift 
distributions, cr p h. 

It remains to be seen how sensitive these conclusions are to 
greater levels of uncertainty in the photometric re dshift distri- 
bution calibration. For example lZhang et aU d2010h consider the 
use of nn and en terms for self-calibrating these uncertainties but 
do not include intrinsic alignments in their main calculations. 
If we assumed e.g. uncertain median redshifts of each individ- 
ual distribution, which is beyond the scope of this paper and 
under investigation elsewhere, we would have obtaine d signifi- 
cantly lower figures o f merit; see lKitching et alJd2008bl) and also 
Bridl e & King! d2007l) who consider this for the case of ee cor- 
relations only. The curve for all correlations in Fig.[5]may hence 
approach the one for ee correlations, especially for <x p h < 0.06. 
Consequently, redshift distributions with a small spread may turn 
out to be even more desirable when taking galaxy number den- 
sity information into account. 
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Fig. 6. The difference o?fom between the FoMtot for a given 
prior value cr pr i or and the fiducial FoMtot, obtained for the least 
stringent, default values of the different priors we apply. Upper 
panel: Effect of tightening the priors on the nuisance parameters 
in the bias terms for ee (dashed line), nn (dotted line), and all 
(solid) correlations. Lower panel: Effect of tightening the priors 
on the slopes of the luminosity function ar", and for the prior 
on the uncertainty of cr ph . Grey curves correspond to priors on 
a w , black curves to those on cr p h- Except for the ee case the latter 
remain very close to zero. As above, ee correlations are shown as 
dashed lines, nn correlations as dotted lines, and all correlations 
as solid lines. Note that ee correlations do not depend on the a w . 
For all curves, the remaining priors are each set to their default 
values of 50 for the bias term parameters, and 10 for priors on 
both or® and cr ph . 



4.3. Dependence on nuisance parameter priors 

On all our nuisance parameters we can expect to have prior in- 
formation to a certain extent, at least by the time large space- 
based cosmological surveys will be undertaken. Here we inves- 
tigate the dependence of the resulting FoMtot on tightening the 
priors on the different sets of nuisance parameters, employing 
our most general configuration with Nk — Nz = 7. We show 
4oM = FoM T oT(o-prior) ~ FoM T ot(ct™* ), i.e. the FoMtot, re- 
ferred to its value for the widest priors we apply as default, as a 
function of the Gaussian prior width <x pl - iol - in Fig. [6] 

For reasons of numerical stability, we have imposed wide 
Gaussian priors on all nuisance parameters, i.e. 10 on both the 
a w and <x ph , and 50 on the parameters within the intrinsic align- 
ment and galaxy bias terms. As the curves clearly indicate, these 
priors are non-informative since all of the curves have flattened 
off by this value. The strongest effect is achieved by tightening 
the priors on the bias terms (note the different scaling of the two 
panels). Prior knowledge on the galaxy bias tremendously im- 
proves constraints by nn correlations, raising the FoMtot close 
to the values of the lensing-only case if the uncertainty of the 



bias term nuisance parameters is of the order 10 . Since in a 
situation with such excellent prior knowledge about galaxy bias, 
the cuts in angular frequency could be much less stringent, the 
increase in FoM would be even more pronounced. 

The FoMtot for ee correlations changes from its value at 
the default configuration to the lensing-only valu^ (compare to 
Table [4]l in the limited range between cr pl -i OT ~ 1 and cr pr i or ~ 
10~ 2 . Thus coarse prior information about intrinsic alignments 
has little effect on the FoM, but, compared to the galaxy bias, a 
model with comparatively moderate precision can already bring 
the FoMtot back close to its optimum. If all available correla- 
tions are used, the internal calibration works well to constrain 
the nuisance parameters, so that the priors only mildly improve 
the FoMtot- However, bear in mind that in spite of our flexible 
parametrisation, the choice of a fiducial base model for both in- 
trinsic alignments and galaxy bias influences the results. Hence 
the base models should be as realistic as possible. It is ques- 
tionable whether our knowledge about the bias terms will ever 
suffice to impose priors on the bias term nuisance parameters 
beyond the approximately correct form of the base models. 

As already discussed in the foregoing section, the global pa- 
rameter crph is excellently constrained by the nn correlations. 
Therefore even tight priors do not have any effect on the FoMtot 
if nn correlations form part of the data vector. Priors below 
fprior = 10~ 2 increase the information content in ee correlations 
marginally. As far as the uncertainty in the redshift distributions 
is concerned, our investigation is still to be regarded as idealistic. 
We defer the joint analysis of shape and number density corre- 
lations in presence of unknown parameters in each individual 
redshift distribution to future work. 

The lack of knowledge in the slopes of the luminosity func- 
tion a® entering the lensing magnification signal is negligible 
in comparison with the effect of the galaxy bias nuisance pa- 
rameters. If all correlations are considered, the a® are well con- 
strained since the improvement in FoM due to priors is below 
the 1 % level. 

4.4. Information content in the individual signals 

The question arises of which signals contributing to the observ- 
able power spectra contain most of the information about cos- 
mology or suffer most strong ly from the uncertai nty in intrinsic 
alignments and galaxy bias. Brid le & Kind (l2007h have studied 
the effect of the two intrinsic alignment terms on ee correlations 
with a parametrisation very similar to ours. In the following we 
consider nn correlations by repeating the Fisher matrix analysis 
twice: once assuming that there is no lensing magnification sig- 
nal, i.e. only gg contributes to nn; and once assuming that there 
is no intrinsic galaxy clustering, i.e. only mm contributes to nn. 

For this calculation we consider constraints from galaxy 
clustering information alone (nn), and assume that galaxy shape 
information (ee and ne) is not used. We show the resulting 
FoMtot for the most flexible model (Nk = Nz = 7) and a model 
with somewhat less flexibility (Nk = Nz = 3) in Fig. [7] On com- 
paring the second, third and fourth columns from the left we see 
that the model with less freedom (solid grey lines) has approx- 
imately the same FoMtot for the complete nn signal as when 
the nn signal is made up of either one of the pure gg and mm 
terms alone. We now examine how this conclusion is changed 



6 In fact we expect the FoM TO t for <x pl .j ol . — > to be slightly above the 
lensing-only case since the intrinsic alignment signals also constrain 
cosmo logy if their form is perfectly known, see also Bridle & King 
(2007). 



15 



B. Joachimi and S. L. Bridle: Simultaneous measurement of cosmic shear and galaxy number density correlations 





60 








50 






H 
O 
, H 


40 




- 




30 








20 




— N K = N Z = 2, 

— N K = N Z = 7 - 



es nn gg mm gg mm 
full set w. Z— cuts w/ o Z— cuts 



Fig. 7. The FoMtot for different subsets of correlations that form 
the nn signal. Throughout, we used N z <om = 10 and cr ph = 0.05. 
The marginalisation has been performed for Nk — Nz = 7 nui- 
sance parameters (black dashed bars) or Nk = Nz = 3 nuisance 
parameters (grey solid bars). The two leftmost columns show the 
full set of ee and nn correlations for reference. The two centre 
columns stand for gg and mm correlations only, employing the 
cuts in angular frequency. The rightmost columns show again gg 
and mm, but here the full information up to ( max = 3000 is con- 
sidered. Note that mm correlations depend neither on intrinsic 
alignments nor on galaxy bias and hence are independent of Nk 
and Nz- 



when increased flexibility is allowed in the galaxy bias model 
(dashed lines). Lensing magnification does not depend on any of 
our bias terms, so that its FoMtot value is the same for the flexi- 
ble model (the dashed and solid lines are on top of each other in 
the fourth column from the left). In contrast, constraints on cos- 
mology from the intrinsic galaxy clustering information alone 
(gg) weaken significantly as greater flexibility is included in the 
galaxy bias model. Specifically, the constraint from the gg signal 
drops by more than 20 in the logarithmic FoMtot when increas- 
ing the number of nuisance parameters per bias term by a factor 
of five (compare the solid and dashed lines in the third column 
from the left). 

The nn correlations, containing both clustering and magni- 
fication, unsurprisingly have a FoMtot in between (second col- 
umn from the left). Hence, if the galaxy bias is well known, the 
galaxy clustering signal dominates the information from the nn 
correlations, but for very flexible bias models lensing magnifica- 
tion does allow more information on cosmology to be extracted, 
as compared to what might be expected if magnification did not 
take place (compare the dashed lines in the second and third 
columns). 

The two right hand columns of Fig. [7] show the FoMtot 
obtained when the full range of angular scales is used in 
the power spectra, thus including information from non-linear 
scales. As already seen in Fig. [2] the clustering constraints im- 
prove strongly when adding the signal from the non-linear 
regime. The magnification contribution remains trustworthy far 
into the non-linear regime, but is of course also affected by the 
/"-cuts. When using the full range of angular scales the FoMtot 
of the mm signal is increased by about 7. The cosmological con- 
straints from intrinsic galaxy clustering (gg) alone are now much 
tighter than from the magnification effect alone (mm), when the 
smaller number of bias parameters are used (compare solid lines 
in the two right hand columns). However, the cosmological con- 
straints from the maximally flexible model are still tighter from 



magnification alone than from galaxy clustering (dashed lines in 
the right hand two columns). 

Note that the constraints from the isolated mm and gg sig- 
nals are hypothetical - the full nn correlations including both 
terms are the only true observables. It is only possible to sep- 
arate the contributions in an approximate fashion when making 
use of their characteristic scaling with redshift, see Fig.[T] But all 
of this available information is already included in our nn results 
when a large enough number of tomographic redshift bins are 
used. 

We now assess the cosmological information available in the 
different subsets of observables when performing the Fisher ma- 
trix analysis for all correlations simultaneously. To this end, we 
split up the summation of the Fisher matrix d42l into three parts, 
corresponding to ee, nn, and ne correlations. We also consider 
pairs of observables e.g. ee, and nn. The Fisher matrices of each 
part are then inverted separately to yield individual parameter 
errors and figures of merit. Thereby we split the total informa- 
tion into subsets of the data vector, which could in principle be 
observed independently. 

However, the figures of merit we compute do not correspond 
to those which one would obtain for an independent analysis of 
ee, nn, or ne correlations because we extract the relevant rows 
and columns of the inverse of the full covariance matrix (l39b 
to insert into the Fisher matrix calculation ( 1421 . Due to the in- 
version, the covariance terms of the different subsets mix, thus 
accounting for the cross-correlations between the observables in 
the different subsets. This is desirable for our purposes because 
these terms add together in the full calculation. Consequently, it 
is possible that the FoM of a subset is larger than the one ob- 
tained for the complete data vector if there are anti-correlations 
with observables of other subsets which are not taken into ac- 
count due to the splitting. Formally speaking, this means that due 
to anti-correlations of different observables, off-diagonal terms 
of the Fisher matrix can have negative entries, which produce 
negative terms in the sum in d42l) . We indeed observe this be- 
haviour for one of the subsets. 

The FoMtot for the subsets of ee, nn, and ne correlations, as 
well as all possible combinations thereof, are given in the upper 
panel of Fig.[8] again for the two parametrisations Nk = Nz = 
{3; 7). As part of the full set of observables, the ee correlations, 
governed by the cosmic shear signal, yield the highest FoMtot 
and contribute most to the cosmological information in the full 
set (second column). They are followed by cross-correlations 
between number density and ellipticity, which have a moder- 
ate FoMtot (fourth column). This occurs in spite of the smaller 
range in angular frequency used for this observable. In addition, 
the least amount of information is lost when switching to the 
more flexible bias model. The subset of nn correlations has by 
far the lowest FoMtot, which becomes even more pronounced 
when the bias model has more nuisance parameters. This hierar- 
chy in FoM is also mirrored in the results for the different com- 
binations of two subsets. 

To get a more explicit measure of the effect of the marginal- 
isation over the galaxy bias and intrinsic alignments nuisance 
parameters, we compute the quantity 

AFoM = In [detf cosm .] -FoMtot (48) 
= - In [det (Feo^.)" 1 ] + In [det (^'^J 

for each of the subsets. The first term in d48b is of a simi- 
lar form as the FoMtot, but the cosmological elements are ex- 
tracted before the inverse or determinant is taken i.e. it does not 
include the marginalisation over nuisance parameters. Hence, 
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Fig. 8. Upper panel: The FoMtot for different subsets of cor- 
relations that are observable from a galaxy survey with galaxy 
shape and number density measurements. The Fisher matrix has 
been computed for the total data set containing all correlations, 
i.e. these results represent the information content of the subsets 
as a part of the total signal. Black dashed bars correspond to a 
marginalisation over Nk — Nz = 7 nuisance parameters, grey 
solid bars to Nk = Nz = 3 nuisance parameters. The labels on 
the abscissa indicate the different combinations of correlations 
used, where en stands for the cross-correlations between num- 
ber density and ellipticity. Lower panel: Difference A FoM, see 
(1481 1. for the same subsets as above. This difference can be under- 
stood as a measure of the depletion of information in the subsets 
due to the marginalisation over nuisance parameters. 



AFoM quantifies the depletion in FoM due to the marginalisa- 
tion. As is evident from the lower panel in Fig. [8] the suscepti- 
bility of the subsets to the nuisance parameters is closely related 
to their contribution to the FoMtot of the full set. Accordingly, 
AFoM is largest for the nn subset while the ellipticity-number 
density cross-correlations have substantially smaller AFoM al- 
though they contain all four bias terms. The effect of intrinsic 
alignments on the ee subset is relatively small when compared 
to the effect of galaxy bias on the nn correlations. 

Since the ee correlations only depend on intrinsic alignments 
(bi and rf) whereas nn correlations only feature galaxy bias (b g 
and b g ), combining them helps little in reducing AFoM, see the 
column labelled 'ee + nn\ Both lines in this column remain 
higher than for the other combinations 'ee + en'and 'nn + en'. 
The ellipticity-number density cross-correlations instead have 
great potential in breaking degeneracies between cosmological 
and nuisance parameters. In particular, adding their subsets to 
the ee correlations further decreases AFoM below the ee-only 
value for the most flexible bias term model. This synergy is pre- 
sumably related to the inter nal calibration between the IG and 
gl signals, as investigated by Zhang] (l2010h . Finally, the full set 



of shape and number density observables clearly calibrates the 
nuisance parameters best as it produces the smallest AFoM for 
both the most flexible and the more rigid bias term models. 



5. Conclusions 

In this work we studied the joint analysis of galaxy number den- 
sity and shape correlations to constrain cosmological parame- 
ters in presence of contaminations by the intrinsic alignment of 
galaxies and the galaxy bias. We considered the signals due to 
gravitational shear, intrinsic shear, intrinsic galaxy clustering, 
and lensing magnification, explicitly computing all possible two- 
point correlations thereof. We introduced a two-dimensional grid 
parametrisation to account for the unknown scale and redshift 
dependence of both intrinsic alignments and galaxy bias. Further 
nuisance parameters were used to describe the uncertainty in the 
width of the photo-z bin-wise redshift distributions and in the 
slope of the galaxy luminosity functions within each bin. 

Our Fisher matrix analysis demonstrates that the simultane- 
ous use of ellipticity correlations, number-density correlations, 
and in particular ellipticity-number-density cross-correlations 
allows for a substantial amount of internal calibration of the bias 
terms. With flexible models that contain in total more than 200 
nuisance parameters we can recover the volume of the error el- 
lipsoid in parameter space when compared to assuming a pure 
gravitational lensing signal and using just the ellipticity correla- 
tion information. The dark energy parameters wq and w a suffer 
more than other parameters on marginalisation over nuisance pa- 
rameters, so that 56 % of the FoMdetf are retained for a Euclid- 
like survey in this most flexible setup considered. The FoMdetf 
for the combined set of shape and number density correlations is 
close to the pure lensing FoMdetf if we choose a model which 
uses about 100 nuisance parameters to describe intrinsic align- 
ments and galaxy bias. Our approach also proves beneficial for 
upcoming ground-based surveys with DES-like survey charac- 
teristics. The benefit is greatest for the more ambitious survey. In 
addition, we assumed the slopes of the galaxy luminosity func- 
tions in each photo-z bin to be unknown nuisance parameters 
and found that they are well calibrated internally from the data. 

The information which we added on top of the standard cos- 
mic shear analysis comes without any extra cost since galaxy 
number density measurements are directly available from imag- 
ing data. Given our encouraging findings, we hence suggest that 
the joint consideration of galaxy shape and number density infor- 
mation could become the standard technique whenever intrinsic 
alignments are suspected to make a significant contribution to 
the cosmic shear signal. 

When interpreting our results, one has to keep in mind that 
our grid parametrisation of the bias terms has limited flexibil- 
ity. On increasing the number of nuisance parameters in the 
grids, the curves for the figures of merit flatten off. However 
we are at present not able to say whether they approach a lower 
limit or continue to decrease for a large number of grid nodes. 
This issue, which is of considerable theoretical and practical 
interest, is currently under investigation, as well as a compari- 
son with the performance of remov al techniques such as nulling 
(iJoachimi & Schneidedl2008ll2009l) . 

Due to the finite number of grid points the resulting angular 
power spectra will not span perfectly the range of possible phys- 
ical models. Moreover, we found that of order 10 nuisance pa- 
rameters in k per grid node in redshift do not suffice to represent 
all relevant functional forms of the galaxy bias. Therefore we 
removed all observables from the parameter estimation which 
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have a significant contribution from the galaxy clustering sig- 
nal at large k where the signal is strongest and the bias least 
known. However, since in the near future we expect to have at 
least coarse, but reliable knowledge about the functional forms 
of the galaxy bias on linear scales and the intrinsic alignment sig- 
nals, the models used in our approach should still yield realistic 
results. 

While the log-linear grid parametrisation is fairly general 
and intuitive, it may not be the most efficient way to represent 
freedom in intrinsic alignments and the galaxy bias. As all bias 
terms originate from physical processes, they are smooth and 
do not feature strong oscillations or isolated peaks. Thus, we 
presume that the bias terms can efficiently be parametrised in 
terms of complete sets of smooth functions such as the Fourier, 
Chebyshev, or Legendre series. Truncating the higher orders of 
these series will only limit the model to represent highly os- 
cillatory or small-scale features. These parametrisations might 
therefore be more effici ent in comprising t he full set of realis- 
tic bias terms (see also Kitc hing et alJl2009h for a given number 
of nuisance parameters, this number in turn being dictated by the 
available computational power. The latter will play a critical role 
when the approach suggested here is performed in a full like- 
lihood analysis with several hundreds of nuisance parameters. 
We add that it might be necessary to furthermore consider the 
bias terms individually for different galaxy types and luminosi- 
ties because it is known that both the intrinsic alignments and the 
galaxy bias vary conside rably with galaxy type and luminosity 
(for recent examples se e Mandelbaum et al. 2006; Hirata et al. 
20071; iMcCracken et al.ll2008HSwanson et al.ll2008t lSimon et al. 
2009tlCresswell & Perciva]||2009HWang et alj|2007l) . 

Redshift information is crucial for discriminating the differ- 
ent signals that contribute to the observables. We investigated 
the dependence of the parameter constraints from the joint set of 
correlations on characteristics of the redshift distributions. We 
confirm the observation by foregoing works that the number of 
photo-z bins needed to retrieve the bulk of information about 
cosmology increases by at least a factor of two when using only 
ellipticity correlations and marginalising over the intrinsic align- 
ment signals. Using the complete set of correlations, the figures 
of merit do not level off so quickly, but continue to increase as 
the number of tomographic bins is increased. This can be ex- 
plained by the narrower kernel in the redshift integrations that 
link the projected and the three-dimensional power spectra in- 
volving galaxy number-density signals. 

Moreover the figures of merit decrease substantially more 
slowly as a function of the photometric redshift dispersion <x ph 
when adding galaxy number density information to the data. 
However, in our marginalisation we only used a single param- 
eter accounting for the uncertainty in the spread of all bin-wise 
redshift distributions. This way the redshift cross-terms of the 
number-density correlations efficiently calibrate the shape of the 
redshift distributions. Therefore this result alone cannot be in- 
terpreted as a potential relaxation of the requirements for photo- 
metric redshift accuracy in cosmological surveys featuring cos- 
mic shear. Conversely, the continuing increase in the figure of 
merit as a function of 7V z bi n makes an even larger number of well- 
separated and compact redshift distributions desirable in the case 
of the joint data set. 

The recent work by |Zhang et al] (1201 Oh uses density- 
ellipticity correlations to self-calibrate photometric redshift pa- 
rameters but a full approach may need to simultaneously deal 
with intrinsic alignments and photometric redshift properties if 
sufficient spectra cannot be obtained to calibrate the photometric 
redshifts independently. In a forthcoming paper we will inves- 



tigate into a more realistic approach including uncertainty and 
outliers in the redshift distributions of each individual photo-z 
bin, as well as the benefits of spectroscopic redshift information 
for a subsample of the galaxy catalogue. 

We emphasise that in the approach suggested here all con- 
sidered correlations help in constraining cosmological parame- 
ters. Hence, none of the contributions is regarded as a systematic 
signal, and the different levels of uncertainty concerning the ex- 
act form of the signals are represented by nuisance parameters. 
By means of the joint analysis we increase the statistical power 
via the cross-calibration abilities of the signals, and reduce the 
risk of undetected systematic effects one faces when consider- 
ing signals such as cosmic shear, galaxy clustering, or lensing 
magnification individually. 

This integrative ansatz is not limited to two-point corre- 
lations, but can be generalised to the three-point level in a 
straightforward way. Future surveys will provide excellent data 
for studying three-point correlations whose exploitation can 
break para meter degeneracies an d improve constraints consid- 
erably (e.g. iTakada & Ja in 2004). Our knowledge about intrin- 
sic alignments and galaxy bias at the three-point level is cur- 
ren tly even more limited than in the two-point case (however 
see ISemboloni ef al. 2008), so that a joint investigation of shape 
and number density observables, including a general parametri- 
sation of the various bias terms, may be an appropriate way 
forward. However, limitations due to computational power will 
most likely play a dominant role in this case. 

One can also think of incorporating further sets of observ- 
ables into the analysis. For instance, planned surveys like Euclid 
will also include a spectroscopic survey of a subset of the galax- 
ies in order to determine the underlying redshift distributions of 
the photo-z bins. The spectroscopic data can be used to mea- 
sure baryonic acoustic oscillations (note that the acoustic peaks 
are not included in our galaxy number-density correlations) and 
galaxy peculiar velocities. The latter allow for the measurement 
of redshift space distortions. Similar to the case considered here, 
the joint analy sis of galaxy num ber density, shape, and velocity 
information (Guzi k et al.l 120101) will efficiently cross-calibrate 
nuisance terms such as the galaxy bias and tighten constraints 
on the cosmological model. 
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Appendix A: Projected lensing magnification signal 

This appendix details the derivation of the enhancement or de- 
pletion of projected galaxy counts njf(fl) due to lensing magni- 
fication. The number density of galaxies n, counted above a flux 
threshold S at angular position and comoving distance x is 
altered by gravitational lensing according to 



n(>S,f K (x)0,x) = 



1 



m,x) 



«0 > 



S 



Jk(xW,x 



(A.1) 



where no is the original galaxy number density, and wh ere yU 
denotes the magnification (Bartelmann & Schneider 2001). One 
assumes that the galaxy luminosity function close to the flux 
limit of the survey can locally be written as a power law, n(> 
S , /k(x)0,x) x S~ a(x \ The slope a depends on the line-of-sight 
distance, or equivalently, redshift, but should not depend on an- 
gular dimensions due to isotropy. However, it is a function of 
the magnitude limit in the observed filter, in this work denoted 
by Him- This dependence is dealt with in Sect. 13.21 but for ease 
of notation we drop r\\ m as an argument of a for the remainder 
of this section. Plugging the power-law form of the luminosity 
function into dA.U yields 



n(>S,f K ( X )0,x) 



■mx) 



(A.2) 



no(>S,f K (x)0,x) 

Again following Bartel mann & Schneideii d2001l) . one can 
approximate the magnification in the weak lensing regime as 
fi « 1 + 2kg- Since kg <s 1, we can in addition do a Taylor ap- 
proximation to arrive at 

n{>S,Mx)0,x) 



n O (>S,f K ( X )0,x) 



{\+2k g {0,x)T {x) - 1 
\+2{a(x)-\)K G {e,x)- 



(A.3) 
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3 - - 



Defining the excess galaxy density contrast due to magnification 
effects as 



n(>S,f K (x)0,x) 



1 



n (>s,f K (x)e, x ) 

= 2 (a(x) - 1) k g (0,x) , 



(A.4) 



one obtains for the corresponding projected density contrast 



f 

Jo 

f 

Jo 



(A.5) 



dXP {i) (x)2(a(x)-l)K G (0, X ) 



a 2 
l 





20 mag 

21 mag 

22 mag 
24 mag 
26 mag 
28 mag 




In exact analogy to the standard derivation of (0 one can 
now insert the relation between the convergence and the three- 
dimensional matter density contrast, 



0.0 0.2 0.4 0.6 

z 



0.8 1.0 



kg(0,x) = 



3// 2 Q n 
~2c r 



V ^OJ^zJ^l WW (A6) 

" s ACr) »<\'> 



which, after swapping the order of integration and the names of 
the integration variables, yields 



f 

Jo 



4^5®Cf)*C/kCf)^) 



Here, we have defined the weight 
3// Q m Mx) 



q {i) (x) 



2 c 2 a(x) 



rXhor 

dx' P ( 'Xx') 
Jr 



(A.7) 



(A.8) 



x fr<y f 2{a( X ')-l). 
fit(x) 

Given that the slope of the luminosity function should be a 
smooth function of comoving distance, a(x) varies only weakly 
over the range of the integration in dA.81 >. being determined by 
the distribution p^(x), which has relatively compact support. 
Hence, the mean value theorem constitutes a good approxima- 
tion, so that we can write 



q®(x)*2(a«>-l)qv>(x) 



(0 



(A.9) 



where we define a (l) to be the slope of the luminosity func- 
tion, evaluated at the median redshift of the photometric bin 
i. Inserting ( |A.9t into dA.7t results in (fT3l ). which we employ 
throughout this work. 

Appendix B: Limber equations of cosmological 
signals 

For reference, we have collected here the Fourier-space Limber 
equations of all projected power spectra that contribute to the 
cosmological signals considered in this work: 
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Fig. C.l. Top panel: Slope of the galaxy luminosity function a 
as a function of redshift, shown for different magnitude limits 
r\i m as indicated in the legend. Bottom panel: Fit parameters a\ 
(solid), «2 (dotted), and 03 (dashed) as a function of ry im . These 
parameters are obtained from the polynomial fit in (1C.21 >. 
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The terminology of both projected and three-dimensional power 
spectra is summarised in Table Q] The weights that enter the 
foregoing equations are the probability distribution of galaxies 
with comoving distance p ( 'Xx)> an d tne ones defined in ( [Tol l, and 
SAM. 



Appendix C: The slope of the galaxy luminosity 
function 

To compute the lensing magnification signal, we need to model 
the power-law slope of the cumulative galaxy luminosity func- 
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Table C.l. Fit parameters for the slope of the luminosity func- 
tion as a function of limiting magnitude ri; m and redshift, see 
(IC2l l and E2). 



j 




bij 


hj 


1 


0.44827 








2 


-1 


+ 1 


+1 


3 


0.05617 


0.19658 


0.18107 


4 


0.07704 


3.31359 


3.05213 


5 


-11.3768 


-2.5028 


-2.5027 



tion at the magnitude limit of the galaxy number densi ty samples 
under consideration. We base our modelling on Blake & Bridle 
(2005) who have determined galaxy redshift distributions for a 
given magnitude lim it, using COMBO- 17 luminosity functions 
for the SDSS r filter (IWolf et al.ll2003h . 

They fitted these distributions with functions of the form (1291 
with /3 = 1.5, using two free parameters, the redshift scaling z 
and the normalisation given by the galaxy surface density So- 
We set Xo and z as a func t ion of survey depth making use of 
Table 1 of Blake & Bridli J2001 and fit a power law to each 
quantity as a function of the limiting magnitude rn m , 

Eo = MM* (C1) 

Z = Z c + Z m (Him - 24) 

where we find good fits using Eo,c = 9.83, rjz = 19, z c = 0.39 
and z m = 0.055. This allows us to extrapolate beyond the range 
of their Table, which stops at nim = 24. 

Note that our definition of a is not the same as the exponent 
of the Sch echter function (sometimes also denoted a, see e.g. 
IWolf et alj |2003. eq. 4). Our a is the negative of the slope of the 
cumulative luminosity function. Therefore for comparison to the 
Schechter function, in the faint galaxy (power-law) regime, one 
must take the negative of our alpha and subtract unity. Note that 
our use of the cumulative luminosity functi on is consistent wit h 
the literature on lensing magnification (e.g. Sch midt et al.l2 009). 
Our typical a values of around 0.5 are also therefore consi stent 
with the luminosity function literature (e.g. lLiu et al .1 120081 find 
Schechter function exponent values between - 1 and -2 depend- 
ing on the spectral type from COSMOS). 

We are interested in the slope of the luminosity function 
<x(z, Him) at the cosmic shear survey magnitude limit. This slope 
is a function of redshift and magnitude limit. From d29b and 
( IC.ll ) we have the number of galaxies as a function of redshift 
and magnitude limit. We convert each magnitude limit into a flux 
limit S and set the number of galaxies above the flux limit equal 
to S~ a (see also AppendixlAt. The resulting curves for a(z, rn m ) 
are shown in Fig. IC.ll upper panel. 

For convenience, and to extrapolate the slope a(z, Him) to val- 
ues z > 1, we now provide a fitting formula . First we expand the 
slope using a polynomial in redshift, with coefficients that de- 
pend on the limiting magnitude. Then we find an approximate 
equation for these coefficients as a function of limiting magni- 
tude. This results in equations for the slope as a function of red- 
shift and magnitude limit in terms of 15 numbers given in Table 
PI 

We fit the slope of the luminosity function as a function of 
redshift with a second-order polynomial 

a(z, Him) = Ol (Him) + 02 (Him) Z + O30"lim) Z 2 . (C.2) 
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Fig. C.2. Percentage deviation of the simplified fit as defined in 
( IC.2b and ( IC.3I ) from the original fits described by (IC.ll ) and 
shown in the upper panel of Fig. IC.ll The contour levels cor- 
respond to the percentages given on the curves. Negative devi- 
ations are indicated by dashed contours. Note that across most 
of the parameter plane the modulus of the deviation is less than 
1%. 

The polynomial coefficients a,- are functions of the limiting mag- 
nitude, and are shown in Fig. IC.ll lower panel. We find that these 
coefficients are in turn well fit by a function of the form 

fli(nim) =b n + b t2 (b t3 ri im - bu) b ' s , (C.3) 

with parameters b, ■ j given in Table IC.ll We chose not to use b t 2 
as a free parameter for the fit, but set it as bi 2 = + 1 , to determine 
the sign of the term in parentheses. By means of ( IC.2l i and (IC.3b 
we have condensed the dependence of a on redshift and r\ lm into 
the 15 parameters summarised in Table IC.ll 

In Fig. IC.2l we plot the relative accuracy of this set of fit for- 
mulae with respect to a(z, Him) as given in Fig. IC.ll lower panel. 
Over the dominant part of the considered parameter space the 
fit formulae provide an excellent approximation, which deviates 
less than 1 % from the original fits (IC.ll ). Thus one can expect 
that within the framework of this approach ( IC.2I ) and ( IC.3b ex- 
trapolate a(z, rii m ) reasonably well to z > 1. Significantly larger 
deviations can only be found for the brightest limiting magni- 
tudes at redshifts z < 0.1, a region of the parameter plane which 
is irrelevant for a competitive cosmological survey. 

We have used and extrapolated the COMBO 17 luminosity 
functions despite the incompleteness beyond R = 23 and relative 
unreliability in the redshift range 1 .2 < z < 2 due to the lack of 
spectral features in the observing bands used. This could be im- 
proved by using deeper data which have infra-red observations, 
for example COSMOS ground and space data and CFHTLS- 
Wide. We emphasise that this is simply used for our choice of 
fiducial model, and within the modelling it is assumed that the 
luminosity function slopes are an unknown function of redshift 
and are marginalised over. We do not expect the choice of fidu- 
cial model to have a large effect on our analysis. Note that the 
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number of galaxies as a function of redshift is inserted into the 
analysis using the assumptions made by the project teams and is 
not linked to our luminosity function calculations. 

It would be more elegant to simultaneously derive the red- 
shift distribution and luminosity function slope from luminosity 
functions derived from data as a function of redshift. However 
this would make it difficult to compare our results with those in 
the literature from the imaging survey project teams, who have 
already chosen a galaxy redshift distribution. We recommend 
that the project teams also make available luminosity function 
slopes to be used in analyses such as these. 

Appendix D: Modelling the effect of dark energy on 
non-linear growth 

In this appendix we summaris e the approach made i n the pub- 
licly available icosmo codeQ (Refreg ier et al.l 12008) to model 
the non-linear evolution of structure in presence of a dark en- 
ergy equation of state w(z) + — 1, which we have adopted for 
th is work. It is based upon a modification of the halofit routine 
of ISmith et alj d2003l) in which the interpolation between open 
and flat cosmological models is determined by the parameter 
/ = Qa/(1 - £2 m ), where / = corresponds to open universes 
without dark energy (we assume Q m < 1 here), and / = 1 to flat 
ACDM models. 

Now the fact is exploited that certain variable dark energy 
models mimic the expansion history of open CDM universes. 
Using tvo = - 1 /3 and w a = in d28b for a flat ACDM model, it is 
readily seen that the same Hubble parameter is obtained as for an 
open CDM universe with identical Q m but without dark energy, 
where (f28T > plays the role of the curvature term. Motivated by this 
coincidence, the interpolation parameter / in the halofit routine 
is replaced by 



where w(z) — wq + w a z/(l + z), and / = Oa/(1 - O m ) as before. 
If a model does not feature dark energy, /' = / = 0. For a flat 
ACDM model the interpolation now takes place between w(z) = 
-1/3, mimicking the case of an open CDM in the original halofit 
(/' = 0), and the cosmological constant w(z) = -1 (/' = 1). 

The p erformance of this simpli stic ansatz has been tested in 
Fig. 10 of Schrabbacket al. (2010), finding fair agreement with 
the fit formula oflMcDonald et al 1(120061) . Note that the fits to the 
simulations by McDonal d et al.l (|2006), which include a depen- 
dence on wo, are not suitable for direct use in this article due to 
the limited range in cosmological parameters, most notably cr%. 



7 http://www.icosmo.org 




(D.l) 



22 



